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Preface 


The  purpose  of  this  study  was  to  explore  the  capability  of  a  new  method,  through¬ 
fold,  as  a  predictor  of  defensible  error  bounds  on  target  loading.  There  were  two  specific 
goals.  First,  to  develop  the  required  derivation  of  a  bounded  definition  for  target  loading. 
Second,  to  develop  a  method  for  implementing  through-fold  for  verification  of  expected 
performance. 

Users  of  current  unfolding  processes  have  a  need  for  theoretically  based  error 
estimates  of  target  loading  in  the  evaluation  of  nuclear  effects  simulations.  However,  due 
to  an  inherent  ill-posedness,  current  methods  lack  the  capability  for  determining  error 
bounds. 

Guidance  and  support  during  this  study  was  received  from  many  individuals.  First, 
I  am  indebted  to  my  faculty  advisor,  Professor  Kirk  A.  Mathews,  for  his  help  and  extreme 
patience  during  this  study.  His  expertise  in  the  field  and  knowledge  of  Mathematica 
proved  invaluable.  Next  I  would  like  to  thank  Capt.  Russell  Daniel  and  Capt.  Dennis 
Miller.  Their  work  on  unfolding  problems  was  crucial  to  my  understanding  of  the  topic. 
Finally  I  would  like  to  thank  my  roommates  Ron,  Mike,  and  Michelle  for  their  help  and  for 
just  putting  up  with  me  during  the  thesis  quarters. 


Kris  G.  Rongone 
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Abstract 

Various  applications  of  the  Fredholm  integral  equation  appear  in  different  fields  of 
study.  An  application  of  particular  interest  to  the  Air  Force  arises  in  determination  of 
target  loading  from  nuclear  effects  simulations.  Current  techniques  first  unfold  the 
incident  spectrum  then  determine  target  loading;  resulting  spectrum  and  loading  are 
assumed  exact. 

This  study  investigates  the  feasibility  of  a  new  method,  through-fold,  for  directly 
determining  defensible  error  bounds  on  target  loading.  Through-fold  uses  a  priori 
information  to  define  input  data  and  represents  target  response  with  a  linear  combination 
of  instrument  responses  plus  a  remainder  to  derive  a  quadratic  expression  for  exact  target 
loading.  This  study  uses  a  simplified,  linear  version  of  the  quadratic  expression. 

Through-fold  feasibility  is  tested  by  comparing  error  bounds  based  on  three  target 
loading  functions.  The  three  test  cases  include  an  exact  linear  combination  of  instrument 
responses,  the  same  combination  plus  a  positive  remainder,  and  the  same  combination  plus 
a  negative  remainder. 

Total  error  bounds  reduced  from  100%  to  35%  in  cases  #1  and  #2.  In  the  third 
test  case  error  bound  was  reduced  to  48%.  These  results  indicate  that  through-fold  has 
promise  as  a  predictor  of  error  bounds  on  target  loading. 
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FEASIBILITY  EXPLORATION  OF  THROUGH-FOLD  AS  A  PREDICTOR 
FOR  TARGET  LOADING  AND  ASSOCIATED  ERROR  BOUNDS 


I.  Introduction 

Background 

Fredholm  integral  equations  of  the  first  kind  apply  in  many  fields  of  study  including 
acoustics,  optics,  geophysics,  and  aerodynamics  (Daniel,  1988:Chl,  1).  Of  particular 
interest  to  the  Air  Force  is  the  measurement  of  pulsed  radiation  emitted  from  nuclear 
weapons  tests.  These  measurements  can  be  used  to  predict  nuclear  weapons  effects  on 
hardware  by  folding  the  incident  spectrum  with  the  hardware  response  function. 

However,  the  Fredholm  integral  equations  are  first  used  to  determine  the  radiation 
spectrum  incident  on  the  equipment.  The  terminology  in  this  study  originates  from  this 
application. 

The  Fredholm  integral  equation  of  the  first  kind  is: 


b 

Y(e)  =  Js(n)  R(e,J\)dt\ 

a 


where 

Y(e)  =  detector  output  signal 
S(r\)  =  radiation  spectrum 
R(e,r\)  =  detector  response  function 


(1) 
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Currently,  pulsed  radiation  detection  is  accomplished  using  a  series  of  overlapping 
detectors  covering  the  energy  range  of  interest.  Thus,  e,  is  discretized  as  i,  the  detector 
index.  However,  practical  aspects  of  data  collection  force  several  limitations  on  a 
detection  system.  These  limitations  include  use  of  a  finite  number  of  detectors,  limited 
detector  resolution,  and  errors  due  to  recording,  transmitting  and  calibrating.  The  exact 
signal  from  each  detector  can  be  represented  with  the  following  simultaneous  equations: 


Yt  =  J  S(E)Ri(E)dE 
0 


where 

Yl  =  the  exact  signal  of  the  ith  detector 
S(E)  =  the  exact  or  actual  spectrum  emitted 
/?,-(£)  =  the  exact  response  of  the  i  detector 


(2) 


Past  efforts  to  resolve  the  incident  spectra  for  determination  of  target  loading 
(nuclear  effects  on  hardware)  have  culminated  in  a  cross-validation  iterative  unfold 
method  (Miller,  1990).  This  method  uses  a  predicted  spectrum  based  on  a  priori  data,  the 
calibrated  response  functions  of  the  detectors,  and  test  data  (i.e.,  measured  detector 
signals)  to  conduct  a  general  deconvolution  of  the  incident  spectrum.  The  unfolded 
spectrum  is  then  smoothed  to  average  out  the  fitting  process  discontinuities.  After 
iterating  the  unfolded  spectrum  to  an  optimal  point  (determined  using  the  cross-validation 
statistical  method)  the  final  unfolded  spectrum  is  used  as  an  approximation  of  the  incident 
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spectrum.  Typically,  the  analyst  folds  this  incident  spectrum  with  the  target  response 
function  to  determine  target  loading.  Unfortunately,  the  unfolded  incident  radiation 
spectrum  is  assumed  to  be  exactly  correct  because  the  iterative  unfold  method  lacks  an 
independent  means  of  determining  how  well  this  unfolded  spectrum  approximates  the 
incident  spectrum. 

Problem 

Since  there  are  an  infinite  number  of  incident  spectra  which  will  produce  the  same 
detector  outputs,  the  traditional  deconvolution  of  the  Fredholm  equations  is  an  ill-posed 
problem.  Coupling  this  ill-posed  problem  with  the  assumption  that  the  resulting  unfolded 
spectrum  is  exact  leads  to  target  loading  estimates  lacking  defensible  error  bounds. 

The  main  purpose  of  this  study  is  to  investigate  the  feasibility  of  a  new  method  for 
determination  of  bounds  (both  upper  and  lower)  on  the  value  for  target  loading.  This  new 
method  (aptly  termed  through-fold)  bypasses  the  intermediate  step  of  determining  the 
spectrum,  and  the  inherent  errors  associated  with  this  step,  to  directly  determine  defensible 
error  bounds  on  target  loading. 

Scope 

This  study  develops  the  through-fold  method  to  estimate  target  loading  due  to  an 
incident  spectrum.  The  through-fold  method  assumes  the  incident  energy  spectrum  is 
predictable.  Since  experimental  data  are  unavailable,  this  study  is  limited  to  simulated 
incident  energy  spectra.  A  set  of  synthetic  test  cases  are  developed  for  an  energy 
spectrum  ranging  from  1  e{}  to  64  E® ,  where  E®  is  an  arbitrary  energy  unit. 
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Using  synthetic  test  cases  for  spectra  and  target  responses,  the  exact  target  loading 
is  calculated  for  each  case  via  a  simple  folding  routine.  Relative  errors  with  respect  to  the 
through-fold  bounds  are  then  compared  to  predict  the  feasibility  of  the  through-fold 
technique. 

Assumptions 

Along  with  assumptions  associated  with  the  through-fold,  additional  assumptions 
are  required  to  be  consistent  with  Miller’s  (1990)  work. 

1.  The  instrumentation  system  consists  of  17  detectors,  6  open  (fluorescer)  response 
functions  and  1 1  closed  (filtered- fluorescer)  response  functions. 

2.  Detector  resolution,  E/AE,  is  1.5. 

3.  The  instrument  response  functions,  RtfE),  are  exact.  The  exact  forms  are  given  in 
Chapter  II. 

4.  Measurement  errors  follow  a  Gaussian  distribution  (Carter,  1989:14-15). 

5.  The  spectrum,  S(E),  is  exact  subject  to  the  estimate  of  uncertainty,  &S(E). 

6.  The  spectrum  has  a  non-negative  flux  constraint. 

7.  The  target  response  function,  T(E ),  has  a  non-negative  constraint. 

Several  other  assumptions  are  required  in  the  derivation  and  simplification  of  the  through¬ 
fold  equations;  their  discussions  are  deferred  until  they  become  relevant. 
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General  Approach 


This  study’s  approach  is  predominantly  analytical.  First,  the  exact  target  loading 
equation  is  developed.  The  equation  derivation  (detailed  in  Appendix  A)  bypasses  the 
intermediate  step  that  determines  the  incident  spectrum,  and  thus  some  of  the  uncertainties 
associated  with  that  process.  Through-fold  directly  determines  a  bounding  expression 
(exact  correction  terms  are  included)  for  the  exact  target  loading.  Upper  and  lower 
bounds  for  target  loading  are  then  expressed  as  linear  programming  problems  which  can 
be  optimized  with  respect  to  assumed  known  bounds  on  the  corrections.  The  elegance  of 
the  through-fold  method  is  in  the  manipulation  of  basic  definitions  for  exact  terms  to 
develop  the  final  expression.  An  example  of  the  definition  of  an  exact  term  is: 

5  =  5  +  85  (3) 

where 

5  represents  an  exact  value, 

5  is  a  predicted  or  measured  value 
85  is  a  correction  term  (error) 

To  incorporate  experimental  data  into  the  equation,  the  exact  target  response  function  is 
represented  by  a  linear  combination  of  the  instrument  response  functions  plus  a  remainder. 
This  is  discussed  further  in  Chapter  II. 

Second,  data  that  would  normally  be  available  to  the  analyst  is  simulated.  These 
data  include  a  predicted  spectrum  and  a  conservative  estimate  of  its  error  (both  based  on  a 
priori  information),  the  detector  response  functions,  and  a  target  response  function.  Once 
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the  input  data  are  developed,  instrument  signals  are  constructed  by  folding  the  simulated 
spectrum  with  the  detector  response  functions.  Measured  signals  are  then  approximated 
by  simulating  measurement  error  with  random  noise  from  a  Gaussian  distribution. 

Third,  the  final  expression  is  optimized  over  the  choice  of  linear  combinations 
available.  This  reduces  the  overall  error  spread  from  what  would  be  calculated  without 
using  experimental  data  in  the  derivation. 

Finally,  the  method  is  tested  to  confirm  that  the  overall  error  spread  can  be 
reduced  via  through-fold. 

Sequence  of  Presentation 

Chapter  II  presents  a  detailed  discussion  of  the  through-fold  theoretical 
development  including  potential  error  sources  and  expectations  for  the  through-fold 
technique.  Chapter  III  covers  implementation  of  the  through-fold  equations  to  include 
input  data,  error  considerations,  and  optimization.  Chapter  IV  discusses  the  through-fold 
test  cases  and  results.  Chapter  V  presents  the  results  of  this  thesis  .  Chapter  VI  presents 
recommendations  for  further  study. 
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II.  Theoretical  Development 


Introduction 

Along  with  introducing  the  basic  concepts  of  the  through-fold  technique,  Section  I 
defined  the  assumptions  and  goals  of  this  study.  This  section  provides  the  theoretical 
background  of  the  through-fold  method.  The  following  topics  are  discussed:  detection 
systems  and  their  associated  response  functions,  sources  of  error,  a  preliminary  derivation 
including  motivation  for  the  approach,  discussion  of  a  more  detailed  derivation,  and 
expectations  of  the  through-fold  technique. 

Detection  Systems 

A  detector  response  function,  RtfE),  is  the  normalized  output  of  the  detector  due 
to  the  portion  of  the  incident  spectrum  at  energy  level  E.  The  response  functions  are 
assumed  to  be  accurately  calibrated  using  an  adjustable  monoenergetic  source. 

In  a  realistic  experiment,  detector  systems  having  multiple  response  functions  are 
used  simultaneously  and  their  results  are  compared  in  an  effort  to  reduce  measurement 
error.  The  detectors  used  to  measure  spectra  of  nuclear  simulations  are  designed  to  have 
an  output  proportional  to  an  incident  spectrum  in  the  energy  range  E  to  E  +AE  given  by 

yi(E)AE  =  Ri(E)S(E)AE  (4) 

where 

yi  (E)AE  =  the  incremental  output  of  detector  i  due  to  the  spectrum  in  the  energy 
range  E  to  E+AE 
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Ri(E)  =  the  response  of  detector  i  to  the  part  of  the  spectrum  at  energy  E . 


S(E)  A E  =  the  amount  of  the  spectrum  that  is  in  the  energy  range  E  to  E  +AE 


As  stated  in  Section  I,  two  different  types  of  detector  systems  having  either  open 
(fluorescer)  or  closed  (filter-fluorescer)  response  functions  are  used  for  this  study.  Both 
detector  response  functions  are  based  on  simplified  versions  (Daniel:  1988,  App-A)  of  the 
detector  response  functions  for  measuring  the  spectra  of  pulsed  x-radiation  presented  by 
Gorbachenko  et  al.  (1976). 

Open  (Fluorescer)  Response  Functions.  The  open  detector  system  uses  fluorescer 
detectors  with  a  zero  response  for  photon  energies  less  than  the  k-edge  of  the  fluorescer. 
Photons  above  the  k-edge  of  the  fluorescer  material  react  with  the  fluorescer,  causing  a 
cascade  of  x-ray  fluorescence.  The  response  decreases  for  energies  above  the  fluorescer 
k-edge  due  to  a  decreased  probability  of  photon  reactions  at  the  higher  energies.  The 
thickness  and  type  of  fluorescer  material  can  be  varied  to  achieve  a  desired  sensitivity. 

The  fluorescer  response  function  used  for  this  study  is 


R°pen(E)  = 


0 


1-exp 


-3 


E  <  Ei° 


E  >  E{ 


where 

R°pen(E)  =  the  open  function  of  the  i^1  detector 
i?f  =  k-edge  of  the  fluorescer  for  the  i^1  detector 


(5) 


8 


E  -  energy 

Figure  1  is  an  example  of  a  typical  open  response  function.  In  this  example,  the  detector 
has  a  k-edge  of  2£1° ,  where  E®  is  defined  as  the  lowest  k-edge  energy  of  the  entire 
detection  system.  All  energy  units  in  this  study  are  given  the  units  of  E®. 


Figure  1.  Typical  Open  Response  Function 


Closed  (Filtered-Fluorescer)  Response  Functions.  The  closed  detector  system 
consists  of  a  filter  and  a  fluorescer  to  shape  an  open  response  between  the  k-edges  of  the 
filter  and  fluorescer  (Carter,  1989:7-9).  Photons  above  the  filter  k-edge  interact  causing 
fluorescence  and  subsequently  degrade  in  energy.  Photons  getting  past  the  filer  and 


having  energies  above  the  fluorescer  k-edge  interact  at  the  fluorescer  in  a  manner  similar 
to  those  of  the  open  detection  system.  Since  the  k-edges  of  the  filter  and  fluorescer  are 
different,  the  shape  of  the  closed  response  system  differs  from  that  of  an  open  response 
system  (Miller,  1990:9-10).  The  filter-fluorescer  response  function  used  for  this  study  is 


nclosed , 

Rt  ( E )  = 


E 

\_ 

E 


1-exp 


-2 


0 


1-exp 


-2 


E 

K  J 


rexP 


exp 


-0.25 


' E C3 


-1.5 


^n3 


E<E- 


E?  <E<  e] 


E>El 


(6) 


where 

Rclosed^E)  _  the  closed  function  of  the  i^  detector 

e]  =  the  k-edge  of  the  filter  for  the  i^1  detector 

Ei  =  k-edge  of  the  fluorescer  for  the  i^1  detector 
E  =  energy 

Figure  2  is  an  example  of  a  typical  closed  response  function.  In  this  example,  the  detector 
has  a  fluorescer  k-edge  of  2 E®  and  a  filter  k-edge  of  4 E® . 
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Figure  2.  Typical  Closed  Response  Function 


Detector  System  Overlap.  For  nuclear  weapons  effects  simulations,  the  array  of 
open  and  closed  detection  channels  is  designed  to  have  overlapping  response  functions  to 
reduce  calibration  and  measurement  errors.  The  k-edges  defining  the  response  regions  of 
the  seventeen  detectors  used  in  this  study  are  listed  in  Table  1.  The  k-edges  are  ratios 
computed  with  respect  to  E® .  The  overlap  of  closed  response  functions  is  shown  in 
Figure  3  (instruments  1-1 1  are  shown  left  to  right).  Figure  4  shows  the  response  functions 
for  the  collection  of  open  detection  systems  (instruments  1-6  are  shown  left  to  right). 

Note  that  the  open  response  detector  systems  completely  overlap  those  of  the  closed 
response  systems. 
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Table  1.  Detector  K-Edge  Ratios 


Detector#  g 


Log-Log  Plot  of  Closed  Response  Functions 


Sensitivity,  0.01 
(arbitrary 
units) 


0.0001 


5  10.  20. 

Energy,  E/E0 


50.  100. 
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Signal  Detection.  Practical  aspects  of  data  collection  force  several  limitations  on  a 
detection  system.  These  limitations  include  use  of  a  finite  number  of  detectors,  limited 
detector  resolution,  and  errors  due  to  recording,  transmitting  and  calibrating.  As 
discussed  above,  pulsed  radiation  detection  is  accomplished  using  a  series  of  overlapping 
detectors  covering  the  energy  range  of  interest.  From  Eq  (4),  the  measured  signal  from 
each  detector  due  to  the  total  spectrum  is  represented  by 


Yi  =  js(E)Ri(E)dE  (7) 

0 

where 

Yi  =  the  ideal  signal  of  the  i^1  detector 
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S  ( E )  =  the  exact  or  actual  spectrum  emitted 
R[(E)  =  the  actual  response  of  the  i^1  detector 

In  practice,  multiple  energy  groups  are  used  to  numerically  approximate  integrals. 
This  technique  sums  over  small  energy  groups  (bins)  —  in  essence,  a  composite  numerical 
quadrature  (Burden  and  Faires,  1993:Ch  4).  Using  the  multi-group  approach  Eq  (7)  is: 

N 

Yi  =  'ZS(Ej)Ri(Ej)AEj  (8) 

7=1 

where 

S  ( Ej )  =  the  exact  spectrum  value  for  the  j^1  energy  group 

A  Ej  =  the  width  of  the  ji^  energy  group 
N  =  the  number  of  energy  groups 

Sources  of  Error 

Two  main  sources  of  error  contribute  to  uncertainties  in  the  spectrum  and  signals: 
measurement  error  (including  transmission,  recording  and  reading  of  the  ideal  signals)  and 
errors  due  to  the  numeric  approximation  of  integrals. 

To  reduce  measurement  error,  overlapping  detector  response  functions  are  used. 

In  addition,  two  detector  systems  are  employed.  The  analyst  accounts  for  these  errors  by 
specifying  a  measurement  standard  deviation  for  each  detector.  These  deviations  are 
based  on  past  experiments  and  detector  calibrations.  Measurement  and  calibration  errors 
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associated  with  the  detection  system  can  be  handled  theoretically  by  applying  normally 
distributed  random  noise  to  simulated  measured  signals  (Carter,  1989:14-15). 

The  second  error  source  is  in  the  mathematical  process,  (i.e.,  the  evaluation  of  the 
integrals  equations  via  numeric  quadrature).  In  this  case,  the  approximation  of  the  integral 
equation  by  a  summation  of  discrete  values,  as  in  Eq  (8),  leads  to  approximation  errors. 
These  errors  can  be  made  negligible  compared  with  other  errors  by  selecting  small  bin 
widths  in  the  summation  process.  Relative  errors  associated  with  the  number  of 
geometrically  spaced  bins  used  in  the  integral  approximation  are  listed  in  Table  2. 


Table  2.  Relative  Error  Due  to  Numeric  Approximation  of  Integral 


15  BINS 

30  BINS 

60  BINS 

120  BINS 

Instrument 

Closed  #1 

0.1193 

0.0046 

0.0011 

0.0003 

Closed  #8 

0.5008 

0.4244 

0.0042 

0.0010 

Open  #1 

0.0066 

0.0016 

0.0004 

0.0001 

Open  #4 

1.0198 

0.0298 

0.0075 

0.0019 

Target  Loading 

0.0040 

0.0010 

0.0003 

0.0001 

Preliminary  Derivation 

The  through-fold  process  allows  calculation  of  bounds  for  the  target  loading 
without  intermediate  calculation  of  the  incident  spectrum.  This  is  accomplished  by 


assuming  exact  target  loading  can  be  represented  by  a  combination  of  known,  predicted 
and  bounded  terms.  Bounds  on  target  loading  are  subsequently  expressed  as  linear 


programming  problems.  The  key  is  to  develop  a  suitable  target  loading  representation. 
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Definition  of  Exact  Values.  Throughout  this  section,  the  terminology  will  remain 
consistent  with  that  presented  earlier.  Summarizing  the  nomenclature: 

tilde,  denotes  an  exact  quantity,  as  opposed  to  measured 

6  =  delta,  denotes  correction  term  determined  from  measurement  error, 

calibration  error,  or  a  conservative  estimate 

th 

Ri(E)  =  instrument  response  as  a  function  of  energy,  i  instrument 

S(E)  =  spectrum  as  a  function  of  energy 

T(E)  =  target  response  as  a  function  of  energy 

Yi  =  instrument  signal,  ith  instrument 

Z  =  target  loading 

NOTE:  The  absence  of  both  tilde  and  delta  indicates  the  quantity  is 

measured  or  calibrated. 

The  derivation  begins  with  a  definition  of  target  loading.  Target  loading  can  be 
viewed  as  the  signal  which  would  have  been  recorded  by  the  target  had  it  actually  been  an 
error  free  instrument.  Target  loading  is  thus  represented  similarly  to  instrument  signal  by: 

oo 

Z  =  J  S(E)  T(E)  dE  (9) 

0 

None  of  the  values  in  Eq  (9)  are  actually  known.  Each  can,  however,  be  approximated  by 
a  measured  or  predicted  value  plus  a  correction  term.  If  the  target  response  function  is 
further  approximated  as  a  linear  combination  of  instrument  responses  plus  a  remainder,  to 
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keep  the  relationship  exact,  then  all  the  terms  on  the  right  hand  side  of  Eq  (9)  can  be 


replaced  with  known,  predicted,  or  boundable  values.  Defining  pertinent  exact  terms, 


S(E)  =  S(E)+5S(E) 

(10a) 

T(E)  =  T(E)+ST(E) 

(10b) 

Ri(E)  =  Ri(E)+mE) 

(10c) 

Yi=Yi+dYi 

(10d) 

Equations  (10)  should  be  read  as  ‘the  exact  value  equals  a  measured  (or  predicted)  value 
plus  a  correction  term’.  The  correction  terms  are  not  known  since  the  exact  values  are 
unknown;  however,  they  can  be  bounded.  The  four  correction  terms  in  Eqs  (10)  can  be 
based  on  conservative  estimates,  calibration  error,  or  measurement.  An  examination  of 
the  approximate  magnitudes  of  these  correction  terms  will  help  simplify  the  derivation. 

Estimation  of  Error  Bounds.  The  exact  correction  terms  significantly  influence  the 
solution.  Although  the  analyst  has  a  priori  information  and  possibly  unfolded  spectra  to 
use  in  his  spectrum  prediction,  the  ill-posed  problem  dictates  a  conservative  bound  be 
placed  on  the  spectrum  correction,  5 S.  For  example,  this  bound  could  be  ±  50%  or  even 
±  100%. 

Measurement  error,  8T,  consists  of  uncertainties  in  the  measured  detector  output. 
This  study  considers  the  cumulative  error  related  to  data  collection  as  measurement  error. 
Although  using  two  types  of  detector  response  functions  helps  minimize  the  measurement 
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error,  its  magnitude  is  still  significant.  This  study  assumes  bounds  on  8F  to  be  ±  15%. 
This  is  consistent  with  previous  work  by  Miller  (1990),  Carter(1989),  and  Daniel  (1988). 

The  calibration  error  associated  with  the  instruments,  8 R,  is  distinct  from 
measurement  error.  Calibration  error  is  an  inherent  uncertainty  in  the  detector  response 
function.  The  experimenter  minimizes  the  magnitude  of  8/?  by  calibrating  the  detection 
system.  This  study  assumes  calibration  error  is  ±  1%  to  ±3%. 

Finally,  there  is  uncertainty  in  the  target  response  function,  8 T.  The  inability  to 
directly  calibrate  the  target  would  appear  to  dictate  a  conservative  estimate  for  8 T. 
However,  this  study  will  only  consider  well  characterized  target  response  functions. 

Accepting  the  above  magnitudes  for  the  correction  terms,  a  simplifying  assumption 
concerning  correction  terms  is  made.  For  the  remainder  of  the  preliminary  derivation,  the 
magnitudes  of  ST  and  8R  are  considered  negligible  compared  to  8Y  and  8S  and  assumed 
equal  to  zero. 

Development  of  Target  Loading.  As  mentioned  in  Chapter  I,  the  key  idea  of  the 
through-fold  technique  is  that  the  target  response  function  can  be  represented  as  a  linear 
combination  of  instrument  response  functions  and  a  remainder.  This  approach  provides 
the  vehicle  for  incorporation  of  test  data  into  the  through-fold  equations.  The  motive 
behind  this  approximation  is  best  illustrated  through  several  examples. 

Motivation.  By  inspection  of  Eqs  (9)  and  (10)  it  becomes  obvious  that 
exact  target  loading  is  driven  by  the  large  spectrum  correction,  SS.  If  the  dependence  on 
8S  can  be  diminished  (via  incorporation  of  test  data)  then  Z  may  have  tighter  bounds. 
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Suppose  the  exact  target  response  is  the  sum  of  two  exact  instrument  responses: 


T(E)  =  R6(E )  +  R7(E)  (11) 

By  Eq  (9),  exact  target  loading  is  equal  to  the  exact  target  response  folded  with  the  exact 
spectrum.  Using  this  definition  and  Eq  (11),  Z  reduces  to  the  combination  of  exact 
instrument  signals  shown  by  the  following  steps: 

Z  =  jS(E)  T(E)dE 
=  \S(E)(r6(E)  +  R7{E))dE 

(12) 

=  JS(£)  R6(E)dE  +  jS(E)  R7(E)dE 
=  Y6  +  Y7 

The  exact  spectrum  has  now  been  replaced  by  the  exact  signal,  which  can  be  approximated 
using  test  data,  measured  signals,  for  Yj  in  Eq  (lOd). 

Next  suppose  that  the  exact  target  loading  can  be  represented  by  a  linear 
combination  of  the  exact  response  of  all  the  instruments, 

f(E)  =  X«i  %£)  <!3) 
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where  the  aj's  are  coefficients  for  the  i  instruments.  The  aj's  can  be  positive,  negative,  or 
zero.  Again,  exact  target  loading  reduces  to  a  linear  combination  of  response  functions: 


Z  =  JS(£)  T(E)dE 
/ 

=  {?(£)  *,(£) 


V  I 


dE 


=  5>  \S{E)Ri{E)  dE 
i 

=!>;  % 


(14) 


Finally,  suppose  that  T  (E)  is  not  quite  a  linear  combination  of  exact  instrument 
responses,  but  has  a  remainder  term  associated  with  it: 

T(E)  =  J^aiRi(E)  +  Ta(E)  (15) 

i 


where  fa  (E)  is  the  remainder  term  parametrically  dependent  on  the  choice  of  coefficients. 
The  expression  for  exact  target  loading  now  becomes: 


Z  =  \S{E)  T(E)dE 


=  jS(E)  £af  *,•(£)  +  TS(E) 


V  i 


dE 


=  2>i  \S(E)  Ri(E)  dE  +  \s(E)  Ta(E)  dE 
i 

=  5>i  Yi  +  jS(E)  Ta(E)  dE 
i 


(16) 
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Although  not  leading  to  an  immediate  solution,  these  examples  provide  insight  into 
the  form  of  the  target  response  function  approximation.  The  benefit  from  this 
approximation  is  the  conversion  of  the  exact  spectrum  to  the  exact  instrument  responses, 
which  allows  for  the  introduction  of  test  data  into  the  through-fold  equations.  One 
obvious  complication  is  the  addition  of  another  unknown  exact  value,  T $  (E) .  However, 
proper  selection  of  the  aj  will  result  in  a  7~  (E)  which  will  tighten  (optimize)  the  overall 
bound  on  Z  .  Conversely,  a  poor  choice  of  a*  will  force  a  value  for  (E)  which  will 
relax  the  overall  bound  onZ  . 

Derivation.  Since  terms  on  the  right  hand  sides  (RHS)  of  Eqs  (10)  are 
considered  known,  they  can  be  substituted  into  the  exact  target  loading  equation, 
eliminating  the  unknowns.  Substituting  Eq  (10b)  into  Eq  (9)  yields 

oo 

Z  =  JS(£)  [T(E)  +  ST(E)]  dE  (17) 

0 


Neglecting  8 T  (for  now): 


Z  =  J  S(E)  T(E )  dE  (18) 

0 

It  now  becomes  beneficial  to  introduce  an  approximation  for  T(E).  Note  that  the 
replacement  is  for  T(E)  and  not  T(E)  as  it  was  in  the  Motivation  section.  Hence  the 
additional  exact  term  from  there  will  not  complicate  the  derivation. 
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Approximating  T(E)  by  a  linear  combination  of  instrument  response  functions  and 
then  making  the  relationship  exact  by  carrying  a  remainder  term  yields  the  following: 


T(£)  =  Xai  +  73(E)  (19) 


where 

th 

a,  =  a  constant  for  the  1  instrument 

Ta  ( E )  =  a  remainder  based  on  the  choice  of  ai,  defined  by  Eq  (19) 


Substituting  the  last  equation  into  the  expression  for  exact  target  loading  and  expanding, 


Z  =  ^aijS(E)Ri(E)dE  +  \s(E)Ta(E)dE  (20) 

i  0  0 


Since  8R  has  been  assumed  zero,  Ri(E)  can  be  replaced  with  (E) .  By  also 
approximating  the  exact  spectrum  with  its  approximation  from  Eq  (10d),Z  becomes 


J  S(E)  Ta  (E)  dE  +  J 5 S(E)  Ta  (E)  dE 
0  0 

(21) 

00 

=  2>z-  f  +  J  S(E)  T5  (E)  dE  +  J 5 S(E)  Ta  (E)  dE 
i  0  0 


The  importance  of  the  derivation  being  performed  in  the  order  shown  by  Eqs  (17)  to  (21) 
is  to  manipulate  Eq  (9)  into  the  correct  form  for  introduction  of  the  exact  instrument 


00 

Z  =^ai\S(E)Ri(E)dE  + 
i  0 
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signals,  Y  .  Experimental  data  (measured  signals)  can  now  be  explicitly  introduced  by 
approximating  the  exact  signal  using  Eq  (lOd) 


Z  =  YJaiYi  +  +  \S{E)Tn{E)dE  +  jdS(E)Ta(E)dE  (22) 

i  i  0  0 

The  final  expression  for  target  loading  is  found  by  writing  the  above  equation 
explicitly  in  terms  of  the  a;  by  replacing  Tg  (E)  with  its  value  from  Eq  (19). 


Z=  +  5>;8^ 

i 

nE)-'£aiRi(E) 


+  J  S(E) 
0 


fSS(E) 


T(E)-^aiRi(E) 


dE 


dE 


(23) 


After  simplifying  and  reversing  the  order  of  integration  and  summation, 

z=  Y,aiYi  +  i*Yi 
i  i 

oo  oo 

-  £ ajs(£)  Ri(E)  dE  -  ^aijbS(E)Ri(E)  dE  (24) 

i  0  *  0 

CO  oo 

+  Js(£)  T( E)  dE  +  J  bS(E)  T{ E)  dE 
0  0 


Having  the  final  analytic  expression  for  the  target  loading,  and  bounds  on  the 
correction  terms,  Eq  (24)  must  now  be  converted  to  a  numeric  one. 
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Conversion  to  a  Numeric  Solution.  Numeric  computational  concerns 


dictate  that  the  integrals  in  Eq  (24)  be  approximated  using  a  summation  of  terms  over  bins 
covering  the  full  range  of  energy.  An  example  of  this  multi-group  approach  was  shown  in 
Eq  (8).  Using  the  approximation  from  Eq  (8),  Z  becomes 


z=  2>,ij  +  Xai  sri 

i  i 

N  N 

-  '£aiJdS(Ej)Ri(Ej)AEj  -  '£ai'£&S(Ej)  R^E j)  AE ]  (25) 

i  7=1  i  7=1 

N  N 

+  ^S(Ej)T(Ej)AEj  +  ^8S(Ej)T(Ej)AEj 
7=1  7=1 


To  make  calculations  easier,  it  is  beneficial  to  define  an  instrument  response  matrix,  R,  a 
spectrum  vector,  s,  an  instrument  signal  vector,  y,  a  target  response  vector,  t,  and  a 
constant  a;  vector,  a,  as 

(R)y  =  Ri (Ej)  AE j 
(s  )j  =S(Ej ) 

(y  )i=Yi  (26) 

(t)j=T(Ej)AEj 
(a)/  =  a,- 

Now  converting  the  exact  target  loading  into  vector  notation  (detailed  in  Appendix  A): 


24 


z 


a(y  +  Sy  -  R  s 
+  s-t  +  8s  t 


(27) 


-  R-8s) 


Equation  (27)  is  an  exact  expression  for  target  loading  because  the  exact  correction 
vectors,  8s  and  Sy,  have  been  carried  throughout  the  derivation.  Unfortunately,  Eq  (27) 
can  not  be  computed  because  the  exact  delta  values  are  unknown. 

Bounds  on  Exact  Target  Loading.  The  value  of  Z  given  by  Eq  (27)  cannot  be 
calculated  because  the  delta  terms  are  not  known  exactly,  only  within  the  bounds  on  the 
correction  terms  discussed  previously.  The  exact  target  loading  can  be  bound  by 


ZL<Z<ZU  (28) 

where  Zu  and  ZL  are  the  upper  and  lower  bounds  respectively  on  target  loading. 

The  upper  (or  lower)  target  loading  bounds  defined  above  can  be  computed  for 
any  fixed  choice  of  a  by  choosing  the  values  for  8s  and  8y  -  both  being  allowed  to  vary 
(within  the  error  bounds)  at  each  energy  bin  —  to  maximize  (or  minimize)  the  target 
loading  expression,  which  is  implicitly  a  function  of  a.  These  bounds  are  computed  by 
solving  the  following  linear  programs  for  the  upper  and  lower  target  loading  bounds 
respectively: 

Zf/(a)  =  Max  [a  (y  +  8y  -  R  s  -  R-8s)+  s-t  +  8s  t] 

8y,  8s 

.  (29) 

Z  (a)  =  Min  [a(y  +  8y  -  R  s  -  R-8s)+  s-t  +  8s  t] 

8y,  8s 
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In  Eq  (29)  8y  and  8s  are  variables  subject  to  the  following  constraints 


Symin  ^Sy^Sy 


max 


^smin  —  Ss  <  8smax 


(30) 


R,  s,  t,  and  y  are  input  data;  and  a  is  the  variable  for  which  the  bound  is  to  be  determined. 
The  goal  functions  of  both  bounds  in  Eq  (29)  are  the  same  and  are  given  by: 


/(8y,8s)  =  a  (y  +  8y  -  R  s  -  R-8s)+  s  t  +  8s  t 


(31) 


The  upper  bound  is  found  by  maximizing  the  goal  function  over  8y  and  8s.  Similarly,  the 
lower  bound  is  found  by  minimizing  the  goal  function  over  the  8y  and  8s.  The  resulting 
values  for  Zu( a)  and  ZL(a)  from  Eq  (29)  represent  the  largest  possible  values  of  Z 
(maximum  and  minimum  respectively)  on  target  loading  for  fixed  set  of  a.  However,  it 
should  be  possible  to  tighten  these  bounds  by  selecting  appropriate  sets  of  a.  The  process 
of  selecting  the  a’s  resulting  in  the  tightest  possible  bounds  is  accomplished  by  optimizing 
the  upper  and  lower  bound  over  all  possible  choices  of  a.  The  result  is  a  classic  minimum- 
maximum  (mini-max)  problem  given  by: 


Zmin  =  Min[zf/(a) 
a 

(32) 

Z^=Max[zL(a)] 

a 
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Z^jn  is  the  minimum  upper  bound  (least  upper  bound)  for  target  loading  and  is  found  by 
minimizing  Zu  (a)  over  a.  Z,^ax  is  the  maximum  lower  bound  for  target  loading,  found 
by  maximizing  ZL  (a)  over  a.  The  outer  loop  of  the  mini-max  problem  is  an 
unconstrained  (the  a  vector  can  take  on  any  values)  optimization  problem. 

Standard  Form.  The  goal  function  in  Eq  (31)  can  not  directly  be  used  in  a  linear 
programming  problem  because  it  is  not  in  standard  form.  The  standard  form  for  linear 
programming  problems  requires  non-negative  variables.  The  bounded  error  terms  are 
therefore  defined  as  the  difference  of  two  positive  numbers 


8  Y  =  8  Y+  -  8  Y 
8S  =  8S+  -  8  S 


(33) 


where  all  the  terms  on  the  RHS  of  Eq  (33)  are  greater  than  or  equal  to  zero. 

The  inner  loops  of  the  mini-max  problems  will  force  one  of  the  terms  on  the  RHS 
of  Eq  (33)  to  zero  and  the  other  to  its  maximum  magnitude  at  each  energy  bin.  The  first 
term  will  be  zero  if  the  worst  case  error  is  negative  while  the  second  term  will  be  zero  if 
the  worst  case  error  is  positive.  Inserting  the  non-negativity  constraints  into  Eq  (31),  the 
goal  function  to  be  maximized  and  minimized  becomes 


/(8y,8s)  =  a 


y  + 


+  s  t 


8y+  -  8y 
+  8s+t  - 


-  R  s 


8s  t 


R-8s+  +  R-8s  1 

J  (34) 


subject  to  the  following  constraints 
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0<Sy+<8ymax 
0<8y“<|-8ymin 
0  -  8s  <8smax 


(35) 


The  goal  function  of  Eq  (34)  remains  linear  in  the  delta  terms  (as  is  the  analytic  solution). 
The  linearity  of  the  goal  function  lends  itself  to  a  number  of  optimization  routines.  It 
should  be  pointed  out  that  this  method  is  much  more  generally  applicable  than  as  used 
here.  As  a  point  of  fact,  the  delta  terms  (8ymax,  8ymim  8smaX,and  8smin)  need  not  be 
constants  (e.g.,  +15%,  -15%,  +50%,  and  -50%),  but  can  vary  in  magnitude  with  the 
energy  bins. 


Differences  In  The  Full  Derivation 

The  preliminary  derivation  assumed  instrument  calibration  errors  and  target 
calibration  errors  were  negligible.  However,  they  are  not  zero  and  must  be  included  in  a 
formal  derivation  of  Z  .  Appendix  A  contains  a  formal  derivation  of  exact  target  loading. 
This  section  highlights  the  main  differences  in  preliminary  and  formal  derivations. 

Processes  associated  with  the  formal  derivation  are  similar  to  those  found  in  the 
preliminary  derivation;  therefore,  the  full  analytic  solution  appears  similar  to  the  shorter 
exact  target  loading  expression  from  Eq  (24).  The  full  analytic  solution  is 
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z=  Y,aiYi  +  Eai5^- 


X  fl/  J  5(£)  /?/  (£)  dE  -  X  J S^(£)  */  (*)  dE 

i  0  i  0 


Xajs(£)  S/?j (£)  d£  -  £aj8S(£)  hRt{E)dE 
i  0  i  0 


(36) 


+  Js(£)  T(E)  dE  +  J8S(£)  T(E)  dE 
0  0 


+  |S(£)8T(£)d£  +  JSS(£)87X£)d£ 

0  0 


One  difference  between  the  two  derivations  is  that  four  new  integrals  containing 
two  additional  error  terms,  8£;(£)  and  §T(E),  appear  in  Eq  (36).  The  presence  of  the 
additional  terms  do  not  in  themselves  change  the  overall  technique  for  finding  a  solution. 
A  more  significant  difference  between  the  two  expressions  is  that  two  of  the  integrals 
now  contain  more  than  one  correction  term.  This  second  difference  results  in  bounds  on 
Z  changing  from  linear  programming  problems  to  quadratic  programming  problems. 

The  non-linearity  of  Z  is  easily  seen  in  the  numerical  approximation  of  Eq  (36): 


Z  = 


fy 


a-i 


+ 


+  8y  -  R  s  -  R-8s 
8R  s  -  8R-8s 
st  +  8st  +  s-St  + 


8s  -St 


(37) 


This  quadratic  equation  is  much  more  difficult  to  solve  than  the  simplified  linear  equation. 
The  remainder  of  this  study  considers  only  linear  dependence  on  error  terms. 
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Through-fold  Expectations 

The  through-fold  method  is  beneficial  only  if  it  produces  tighter  bounds  than  can 
be  calculated  based  on  estimated  error  bounds.  As  discussed  under  Motivation,  how  well 
the  method  works  depends  on  the  choice  of  aj.  A  good  choice  of  aj  will  result  in  a  value 
for  (E)  which  tighten  bounds  on  Z  .  Conversely,  a  poor  choice  of  a;  could  result  in 

large  bounds  on  Z  .  Since  T^(E)  is  a  remnant  of  the  through-fold  method,  it  can  be 
viewed  as  error  (the  correction)  associated  purely  with  the  technique.  How  well  through¬ 
fold  performs  is  therefore  linked  to  how  well  (E)  is  optimized.  Optimization  over  a 

does  not  mean  minimization  of  T^(E) ,  but  in  choosing  T^(E)  to  tighten  bounds  on  Z  . 

Conservative  Estimate.  A  conservative  estimate  on  the  exact  target  loading  error 
bounds  can  be  calculated  prior  to  performing  an  experiment.  This  error  would  be  based 
solely  on  propagation  of  expected  errors.  By  utilizing  the  inner  loop  of  through-folds’ 
mini-max  routines,  the  conservative  value  can  be  calculated.  This  is  done  by  computing 

ZL{ a)  and  Zu  (a)  with  all  a;  set  to  zero,  eliminating  instrument  influence  on  the  solution. 

How  accurately  through-fold  techniques  calculate  the  conservative  value  depends  on  the 
number  of  bins  used  in  the  numerical  integration. 

Introduction  of  Experimental  Data.  The  introduction  of  experimental  data  into  the 
through-fold  routines  will  help  optimize  T^(E)  and  therefore  tighten  bounds  on  the  exact 
target  loading.  The  reason  for  this  is  seen  by  comparison  of  two  equations.  The  first 
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equation  is  found  from  a  naive  approach  to  the  original  exact  target  loading  expression, 

Eq  (9).  Starting  with  Eq  (9)  and  remembering  bT(E)  has  been  defined  as  zero  yields: 

oo 

Z  =  j  S(E )  T(E)  dE 
0 
oo 

=  js(E)  T(E)dE 
0 

(9b) 

oo 

=  J[S(£)  +  8S(£)]  T(E)dE 
0 

oo  oo 

=  J  [S(E)  T(E)dE]  +  J  [bS(E)  T(E)  dE] 

0  0 


The  second  equation  is  Z  prior  to  simplification,  a  variation  of  Eq  (22): 


z=E«iii 

+ 

i 

i 

=  1^ 

+ 

i 

i 

oo  oo 

+  J  S(E )  Ts  ( E )  dE  +  J  5 S(E)  Ta  (E)  dE 
0  0 

(22b) 

oo 

+  \S{E)Ta{E)dE 
0 


As  seen  by  Eqs  (9b)  and  (22b),  the  only  error  associated  with  the  naive  approach  is 
the  a  priori  spectrum  error  (estimated  at  ±  50%).  After  applying  through-fold  techniques 
the  spectrum  error  is  replaced  with  measurement  error  (estimated  at  ±  15%)  propagating 
into  the  choice  of  a  plus  the  error  associated  with  the  through-fold  method,  Ta  (E) , 

propagating  into  the  exact  spectrum.  As  discussed  previously,  the  proper  choice  of  ai  will 


31 


optimize  T^(E) .  If  T$(E)  is  driven  towards  zero,  then  the  resulting  error  is  driven 
towards  the  measurement  error  (assuming  propagation  of  measurement  error  is  not 
significant).  However,  if  Tg{E)  can  not  be  driven  low  enough,  additional  error  associated 

with  the  integral  will  counteract  some,  or  even  all,  of  the  difference  between  hS  and  87. 

It  is  important  to  realize  that  the  optimal  T^(E)  may  not  always  be  the  minimum. 

If  all  aj  are  found  to  be  positive,  or  all  negative,  then  propagation  of  the  measurement 
errors  will  be  such  that  the  optimal  7~  (E)  is  as  small  as  possible.  However,  if  the  sign  of 
the  a;  vary,  then  propagation  of  measurement  errors  may  dictate  that  the  optimal  (E) 
not  be  its  minimum. 

Incrementing  Experimental  Data.  As  the  number  of  instruments  used  in  the 
through-fold  increases,  the  error  bounds  are  expected  to  tighten.  This  occurs  because  the 
additional  test  data  more  accurately  define  the  incident  spectrum  implicitly  contained  in  the 
instrument  measurements.  But  more  importantly,  with  the  addition  of  each  instrument,  the 
routine  is  given  two  additional  degrees  of  freedom  for  determining  the  aj.  The  theoretical 

minimum  error  bound  is  equal  to  the  measurement  error;  however,  it  is  unlikely  that  this 
minimum  will  be  reached. 

There  are  two  primary  reasons  that  the  optimized  error  bound  is  not  expected  to 
reach  the  theoretical  minimum  error  bound.  First,  as  discussed  previously,  all  the  a;  are 

not  necessarily  the  same  sign.  Second,  the  target  response  function  is  not  a  linear 
combination  of  the  instrument  response  functions. 
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Ill  Computer  Implementation 


Introduction 

Chapter  II  provided  the  theoretical  background  for  understanding  the  through¬ 
fold  concept.  This  section  covers  implementation  of  the  through-fold  equations.  Included 
in  this  discussion  are  the  implementation  platform,  input  data,  error  considerations, 
optimization  routines,  and  the  final  computer  code. 

Platform 

The  computer  program  Mathematical  ran  on  a  personal  computer,  was  chosen  for 
through-fold  implementation.  Along  with  modest  built-in  optimization  routines  and  the 
capability  to  handle  the  vector  equations  from  Chapter  n.  Mathematical  notebook 
interface  allows  creation  of  interactive  files  for  future  use. 

Input  Data 

As  discussed  in  Chapter  n,  a  variety  of  input  data  is  required  to  perform  a  linear 
optimization  of  Eq  (34).  The  data  normally  available  to  the  analyst  from  theoretical 
treatment  of  the  physics  and  from  past  tests  will  be  simulated.  This  includes  the  a  priori 
spectrum  and  instrument  response  functions,  the  target  response  function,  and  all 
associated  error.  This  section  will  define  the  data  used  for  the  original  validation  of  the 
method.  Instrument  response  functions  and  required  errors  were  defined  in  Chapter  II. 

Spectrum.  A  continuous  normalized  Planckian  spectrum  was  chosen  to  simulate 
the  incident  spectrum.  The  equation  for  a  one  temperature  Planckian  spectrum  is 
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S(E)  =  sf 


15 


fcT  % 


\ 


m 


where  E  is  energy,  kT  is  the  temperature  of  the  radiation  spectrum  (units  of  fif),  and  sf 
is  a  scaling  factor  for  the  radiation  spectrum.  A  convenient  single  temperature  validation 
spectrum  (kT  =  l^f)  and  scaling  factor  of  one  was  chosen  for  validation.  Figure  5  is  a 
plot  of  the  validation  spectrum. 


Figure  5.  Validation  Planckian  Spectrum 
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Target  Response  Function.  A  Gaussian  distribution  function  was  chosen  to 
approximate  the  target  response  function,  T(E),  for  the  validation  case.  The  equation  for  a 
Gaussian  distribution  is  given  by: 


T(E)  = 


1 

o  -JT% 


(39) 


where  E  is  energy,  (l  is  the  mean  of  the  distribution,  and  a  is  the  standard  deviation  of 
the  distribution.  A  Gaussian  distribution  with  a  mean  of  zero  and  standard  deviation  of 
15  Ei  was  chosen  for  the  validation  target  response  function.  The  choice  of  test  target 
response  functions  will  be  discussed  Chapter  IV.  Figure  6  is  a  plot  of  the  validation  target 
response  function. 
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Instrument  Signal.  Exact  instrument  signals  are  constructed  by  folding  the 
simulated  spectrum  with  the  detector  response  functions.  The  folding  is  accomplished  via 
the  numerical  approximation  of  Eq  (7)  and  by  assuming  that  the  a  priori  spectrum 
definition  and  calibrated  instrument  response  functions  are  exact.  Calculation  of  the  exact 
signal  is  simplified  in  Mathematica  by  using  the  matrix  and  vector  definitions  from  Eq 
(26).  The  exact  signal  can  then  be  calculated  by: 

Y=  R  S  (40) 

All  matrix  and  vector  functions  from  Eq  (26)  are  evaluated  at  the  geometric 
centers  of  geometrically  spaced  energy  bins.  Bin  spacing  and  centers  were  chosen  to 
correspond  to  the  geometric  spacing  of  the  instrument  response  functions.  Geometrically 
spaced  bins  result  in  the  majority  of  bins  (60%)  within  the  first  quarter  of  the  energy 
range.  Geometrically  spaced  energy  bins  are  desirable  for  this  application  because  of  the 
higher  resolution  provided  to  the  lower  k-edge  detection  systems  than  is  possible  using 
linear  bin  spacing  for  the  same  computational  costs  (i.e.,  the  same  number  of  bins).  The 
advantage  is  due  to  all  the  interesting  details  of  the  chosen  incident  spectrum,  target 
response  functions,  and  instrument  response  functions  occurring  at  the  lower  energies. 

The  Mathematica  notebook  allows  for  computation  of  linearly  spaced  bins  with  linear 
centers  for  cases  having  considerable  detail  at  higher  energies.  Candidate  test  cases  for 
linear  spacing  include  high  temperature  Planckian  spectra. 
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Error  Considerations 


Numerical  Approximation.  Geometric  energy  bin  spacing  is  used  to  reduce  errors 
in  approximating  integrals  with  summations.  As  shown  in  Table  2,  use  of  a  120  bin  grid 
produces  relative  errors  of  less  than  10  in  approximating  the  required  integrals.  To  save 

computational  costs,  the  validation  cases  used  30  bins,  producing  a  relative  error  that 
varies  with  instrument  as  shown  in  Table  2.  Unfortunately,  the  relative  errors  associated 
with  the  instruments  increase  from  0.001  to  0.42  as  the  fluorescer  k-edge  increases.  The 
high  relative  error  associated  with  closed  instrument  #8  would  normally  be  unacceptable; 
however,  the  spectrum  value  associated  with  the  operational  energy  range  of  this 
instrument  (E  >  1 1.5  E\ )  is  approximately  10’3.  Although  the  total  signal  generated  by 

this  instrument  is  highly  uncertain,  its  magnitude  is  too  small  (10  compared  to  10  for 

lower  numbered  instruments)  to  significantly  influence  the  through-fold  results. 

Noise  Simulation.  In  calculating  the  exact  instrument  signals  (Eq  39),  the 
calibrated  instrument  response  functions,  Ri(E),  are  assumed  exact.  However,  the  exact 
target  loading  equations  in  Chapter  II  require  measured  signals.  This  study  simulates 
measurement  error  by  the  addition  of  random  noise.  These  errors  are  assumed  normally 
distributed  and  determined  by  randomly  sampling  a  Gaussian  probability  distribution. 

This  error  approximation  was  performed  using  a  Mathematica  statistical  package. 
A  single  random  value,  based  on  the  instrument’s  standard  deviation,  was  selected  from  a 
normal  distribution.  In  keeping  with  previous  unfolding  studies,  the  instrument  standard 
deviation  is  taken  to  be  15%  and  includes  the  totality  of  measurement  related  errors. 
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Optimization  Routines 


Built  in  Mathematica  routines  were  used  for  the  two  step  mini-max  optimization 
discussed  in  Chapter  0.  The  first  routine  requires  maximizing  and  minimizing  capabilities 
for  the  linear  target  loading  function  over  the  bounds  on  the  error  terms.  Mathematica 
provides  two  vehicles  (LinearProgramming,  and  ConstrainedMax  /ConstrainedMin)  for 
computation  of  global  maxima  and  minima  in  linearly  programming  problems.  The 
constrained  commands  were  chosen  because  they  easily  implement  variable  sized  problems 
and  automatically  calculate  the  constrained  solution  based  on  the  optimal  constraints. 

The  second  optimization  routine  requires  minimization  of  all  possible  maximum 
values  and  maximization  of  all  possible  minimum  values  to  reduce  the  error  band  between 
the  two  cases.  Mathematica  provides  a  minimization  routine  for  unconstrained  problems. 
The  command  FindMinimum  was  used  to  optimize  the  solutions.  The  upper  and  lower 
bounds  were  allowed  to  use  different  sets  of  aj  to  optimize  the  bound.  See  Appendix  B 

for  the  complete  Mathematica  notebook  implementing  these  commands. 

Performance  of  optimization  routines  is  often  linked  to  the  initial  values  provided 
as  starting  points  to  the  program.  This  is  a  concern  with  Mathematical  FindMinimum 
command  as  it  is  not  guaranteed  to  return  a  global  minimum  value.  As  used  with  the 
through-fold  equations,  the  FindMinimum  command  requires  the  input  of  the  first  two 
points  from  which  to  perform  a  search  for  a  local  minimum.  A  method  of  steepest  descent 
on  the  surface  (function)  is  used  by  the  routine  to  find  the  minimum. 

In  development  of  through-fold  a  number  of  tests  were  performed  to  evaluate  the 
sensitivity  of  the  method  to  initial  guesses.  Although  not  performed  in  a  systematic 
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manner,  testing  of  initial  values  around  (and  away  from)  known  minimums  was  conducted. 
The  optimization  routines  were  found  to  be  relatively  insensitive  to  variations  in  the 
starting  values. 

ThruFold.ma 

ThruFold.ma  is  the  final  version  of  the  through-folding  code.  A  printed  copy  of 
ThruFold.ma  is  provided  at  Attachment  B.  The  computer  code  is  a  Mathematica  for 
Windows  version  2.2  notebook  which  optimizes  bounds  on  exact  target  loading  (linear 
dependence  on  error  bounds  only).  The  code  is  written  in  three  sections:  1)  selection  of 
input  data,  2)  evaluation  of  vectors  and  matrixes  and  3)  optimization. 

The  first  section  allows  the  user  to  select  from  a  number  of  incident  spectra,  target 
response  functions,  and  error  terms  through  evaluation  of  appropriate  cells.  The  user  is 
then  allowed  to  select  the  number  of  instruments,  number  of  bins,  bin  spacing  (geometric 
or  linear),  energy  range  to  be  covered,  and  the  fluorescer  k-edge  of  the  first  detector. 
Finally  the  instrument  response  functions  are  defined.  The  second  section  of  ThruFold.ma 
defines  vectors  for  the  exact  and  measured  signals,  spectrum,  target  responses,  all  error 
approximations,  and  the  instrument  response  matrix.  The  final  section  of  the  notebook 
defines  the  exact  target  loading  using  the  final  form  of  the  linear  programming  goal 
function,  Eq  (34),  develops  equations  for  maximum  and  minimum  values,  and  provides 
examples  for  calculation  of  the  final  bounds.  Default  values  fit  the  validation  test  case. 
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IV  Validation  and  Testing 


Introduction 

In  Chapter  IE  the  implementation  of  the  through-fold  methodology  was 
discussed.  In  this  Chapter  the  validation  cases  and  test  cases  are  presented  and  compared 
to  the  expectations  from  Chapter  II. 

Validation 

A  validation  case  was  constructed  to  ensure  proper  operation  of  the  through-fold 
routine.  There  are  two  validation  goals.  The  first  validation  goal  is  to  ensure  that 
through-fold  produces  the  expected  conservative  error  when  all  a;  are  set  to  zero.  The 
second  goal  is  to  verify  that  error  bounds  decreased  incrementally  as  additional 
instruments  were  introduced  to  the  through-fold. 

Proper  verification  of  the  first  goal  involves  comparing  the  through-fold  results 
with  the  exact  target  loading.  As  mentioned  in  Chapter  II,  the  accuracy  of  the  results  are 
dependent  on  the  number  of  bins  used  to  approximate  the  integrals.  Comparisons  will  be 
complicated  by  this  difference  if  not  using  a  substantial  number  of  bins.  For  comparison 
purposes  only,  a  pseudo-exact  target  loading  is  assumed  equal  to  the  dot  product  between 
the  spectrum  and  target  loading  vector.  Vector  dimensions  are  based  on  the  number  of 
bins  used  for  a  given  test.  This  assumption  alleviates  numerical  approximation  error  when 
all  are  set  equal  to  zero. 
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Validation  Case.  The  input  data  for  validation  were  defined  in  Chapter  III  to  be  a 


broad  Gaussian  target  response  function  (mean  =  0,  standard  deviation  =  15),  a  one 
temperature  Planckian  spectrum  (temperature  =  1 E®,  scaling  factor  =  1),  and  30 
geometrically  spaced  bins  covering  an  energy  range  from  1 E®  to  64 E® .  Each 
geometrically  spaced  bin  has  a  width  ratio  (upper  edge  to  lower  edge)  equal  to  V2  . 
Measurement  error  has  been  introduced  to  the  exact  signals  necessitating  the  inclusion  of 
them  in  the  constrained  optimization  routines.  The  validation  case  uses  three  test  data 
variations:  1)  only  closed  instrument  signals,  2)  only  open  instrument  data,  and  3)  signals 
from  both  closed  and  open  instruments  in  the  through-fold. 

Validation  Results.  Results  of  the  validation  method  are  encouraging.  All  three 
variants  determined  the  error  computed  when  all  ai  equal  zero  to  be  ±50%.  Table  3  lists 
the  results  (maximum  and  minimum  bounds  on  exact  target  loading)  from  the  first 
validation  variant-only  closed  instrument  signals.  The  instruments  used  in  each  test  are 
chosen  sequentially  based  on  the  lowest  fluorescer  k-edge  (e.g.,  No.  Insts.  =  4  indicates 
the  first  four  closed  instruments  (left  to  right)  as  shown  in  Figure  3). 


Table  3.  Validation  Results  (Bounds  on  Target  Loading),  Closed  Instruments 


No.  Insts. 

0 

2 

4 

6 

8 

10 

11 

ZU. 

0.03696 

0.03514 

0.03971 

0.03457 

0.03452 

0.03452 

Zjmin 

Pseudo  Exact  Z 

0.02464 

0.02464 

0.02464 

0.02464 

0.02464 

0.02464 

0.02464 

z l 

0.01232 

0.01525 

0.01849 

0.02009 

0.02046 

0.02030 

_ 

■^max 

Upper  Bound,  %  + 

50.0 

42.6 

61.2 

40.3 

40.1 

40.1 

— 

Lower  Bound,  %  - 

50.0 

38.1 

25.0 

18.4 

17.0 

17.6 

— 
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The  bounds  on  target  loading  are  shown  as  a  percentage  relative  to  the  pseudo- 
exact  target  loading.  The  bounds  generally  decrease  from  the  naively  computed  maximum 
of  ±50%  with  no  instruments.  It  is  interesting  that  the  upper  and  lower  bounds  are  not 
driven  to  the  same  value.  This  is  expected  as  separate  optimization  of  upper  and  lower 
bounds  allow  each  to  have  its  own  set  of  a*  and  hence  different  optimal  values. 

The  optimization  routines  were  unable  to  compute  bounds  when  all  eleven  closed 
instruments  were  used.  Although  unfortunate,  this  is  not  unexpected  due  to  the  extremely 
small  signals  associated  with  the  last  few  instruments.  Figure  3  indicates  the  last  closed 
instrument  has  a  maximum  response  less  than  0.01  at  an  energy  of  45  iEr®.  The  spectrum 
value  associated  with  this  energy  is  ~10  16 The  Mathematica  optimization  routine  could 
not  handle  the  small  variations  associated  with  this  instrument. 

Table  3  clearly  shows  that  the  optimal  maximum  error  bound  increased  beyond  the 
naive  value  when  computed  using  four  instruments.  This  was  unexpected  as  the  routine 
should  have  performed  no  worse  than  when  using  two  instruments.  As  a  check  of  the 
expected  insensitivity  on  the  starting  values  provided  to  the  FindMinimum  command  (see 
Chapter  3,  Optimization  Routines),  further  testing  was  performed  on  this  anomaly. 

In  all  14  additional  tests  were  run.  The  first  five  tests  varied  the  starting  values  of 
all  four  instruments  while  the  remaining  tests  concentrated  on  the  starting  location  of  only 
one  instrument.  Four  of  the  first  five  tests  resulted  in  relative  upper  bounds  near  61%, 
supporting  the  previous  conclusion,  relative  insensitivity  to  starting  positions.  The  fifth 
test  resulted  in  a  bound  of  32%,  less  than  the  naive  case,  as  expected,  but  also  less  than 
other  closed  instrument  upper  bounds,  this  was  unexpected  and  lead  to  additional  testing. 
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The  second  set  of  nine  tests  varied  only  the  starting  locations  of  the  fourth 
instrument.  The  fourth  instrument  was  chosen  because  it  was  the  only  instrument  for 
which  initial  values  were  changed  when  the  32%  bound  was  achieved.  Eight  of  these  nine 
test  found  a  relative  bound  near  32%.  There  does  not  appear  to  be  a  set  pattern,  but  as 
long  as  one  (1)  was  not  used  as  the  first  of  the  two  initial  starting  locations  (both  0.5  and 
1.5  work),  then  the  minimum  upper  bound  was  less  than  35%.  The  second  of  the  two 
guesses  did  not  appear  to  be  a  driver  for  the  optimization  routine. 

These  results  indicate  that  there  may  be  two  local  minima  for  the  function. 
However,  since  starting  values  above  and  below  one  (1)  can  be  chosen  that  result  in  a 
32%  bound,  the  two  apparent  minimum  values  may  actually  indicate  a  slowly  varying 
minimum.  Results  from  the  additional  testing  indicate  proper  performance  of  the  through¬ 
fold  method  requires  user  artisanship.  However,  the  importance  of  the  user  involvement 
may  only  be  a  factor  due  to  the  limited  capabilities  of  the  minimization  routine  (a  steepest 
decent  method)  utilized  by  the  FindMinimum  command. 

Another  possible  cause  for  the  observed  erratic  behavior  is  inadequate  resolution 
of  the  instrument  response  functions  as  provided  by  the  30  bin  structure.  This  type  of 
inadequacy  could  result  in  some  closed  instruments  becoming  indistinguishable  in  their 
response  from  the  open  instruments  responses.  However,  the  exact  cause  for  this  anomaly 
remains  unexplained  beyond  the  non-robustness  of  the  Mathematica  kernel. 

A  final  observation  from  validation  case  #1  is  that  the  total  relative  bound  reduced 
from  a  naive  value  of  100%  to  57%  (49%  after  additional  testing)  utilizing  only  closed 
instrument  responses. 
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Validation  results  from  the  open  instruments  are  listed  in  Table  4.  The  results  are 
comparable  to  those  for  closed  instruments.  Naive  error  equals  ±  50%,  bounds  decreased 
as  instruments  were  added,  and  the  minimization  routine  failed  for  the  last  instrument. 
Non-robustness  of  the  optimization  routines  are  again  expected  to  be  the  culprit . 


Table  4.  Validation  Results  (Bonds  on  Target  Loading),  Open  Instruments 


No.  Insts. 

0 

1 

2 

4 

5 

6 

zv. 

0.03696 

0.03576 

0.03131 

0.02970 

0.02970 

0.02970 

^min 

Pseudo  Exact  Z 

0.02464 

0.02464 

0.02464 

0.02464 

0.02464 

0.02464 

zL 

0.01232 

0.01332 

0.01408 

0.01639 

0.01639 

_ 

Upper  Bound,  %+ 

50.0 

45.2 

27.1 

20.5 

20.5 

20.5 

Lower  Bound,  %- 

50.0 

46.0 

42.9 

33.5 

33.5 

— 

There  are  two  differences  between  closed  and  open  instruments  results.  First,  the 
open  responses  correspond  to  tighter  overall  bounds  than  for  closed  instruments  (53% 
compared  to  57%).  The  tighter  bounds  are  predominantly  due  to  a  large  reduction  of  the 
upper  bound.  Second,  the  overall  lower  bound  was  optimized  using  closed  instruments 
while  the  upper  bound  was  optimized  using  the  signals  from  open  instruments.  This  last 
difference  is  consistent  with  allowing  separate  optimizations  based  on  instrument  type. 

The  third  validation  case,  use  of  open  and  closed  instrument  responses  in  the 
through-fold  failed.  The  optimization  routines  were  unable  to  converge  to  a  solution. 
Two  possible  reasons  are  inadequate  bin  spacing  and  hence  insufficient  instrument 
resolution  and  the  non-robustness  of  the  optimization  routines.  However,  the  routine 


should  have  been  capable  of  producing  results  no  worse  than  the  best  combination  of  the 
two  previous  cases  (i.e.,  an  upper  bound  of  ±  20.5  and  a  lower  bound  of  ±  17.5%). 

With  the  exception  of  the  anomaly  observed  in  the  4-closed  channel  results,  the 
validation  test  verified  all  expectations.  Even  though  some  limitations  of  the  optimization 
routines  were  found,  through  proper  selection  of  the  ai  vector  and  prudent  choice  of 
instrument  response  systems  used  in  the  through-fold,  the  overall  error  bound  was 
decreased  from  100%  to  38%  by  the  through-fold  technique.  Remember  the  theoretical 
minimum  was  30%  in  the  validation  cases  due  to  the  built-in  ±15%  measurement  error. 

Testing 

A  few  changes  were  made  to  the  optimization  problem  between  the  validation  and 
test  cases.  In  an  effort  to  save  computational  time  the  number  of  bins  was  reduced  to  20. 
In  order  to  counteract  the  subsequent  decrease  in  the  accuracy  of  the  numerical 
calculations,  the  energy  range  was  decreased  from  64  to  40  energy  units.  The  resulting 
numerical  approximation  error  was  verified  to  be  comparable  to  the  original  error  using  30 
bins  (0.005  to  0.04  relative  error  —  with  the  exception  of  high  numbered  closed 
instruments).  This  strategy  worked  in  reducing  the  error  due  to  the  two  factors  discussed 
previously,  geometric  bin  spacing  and  the  low  spectrum  values  beyond  energy  unit  15 E® . 

Only  12  of  the  original  17  detector  systems  are  used,  eight  closed  and  four  open 
detection  systems.  There  are  two  reasons  for  this  change.  First,  by  decreasing  the  energy 
range  some  instrument  response  functions  fall  predominantly  outside  the  area  of  concern. 
Second,  the  validation  (along  with  all  other  trials)  resulted  in  the  optimized  error  bounds 
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being  reached  prior  to  incorporation  of  the  eighth  closed  and  fourth  open  instrument.  It 
was  assumed  that  the  test  cases  would  behave  in  a  similar  manner  as  long  as  the  target 
response  function  did  not  significantly  increase  (compared  to  the  validation  case)  at 

energies  above  20 E®.  Lastly,  based  on  validation  results  from  case  3,  optimization  was 
not  performed  using  signals  from  open  and  closed  detection  systems  simultaneously. 

Test  Cases.  Three  test  cases  were  constructed  to  study  the  behavior  of  through¬ 
fold,  specifically  the  influence  of  T^(E)  on  the  solution.  All  three  cases  are  based  on  a 

target  response  function  that  is  a  linear  combination  of  instrument  response  functions. 
Since  through-fold  is  designed  to  optimize  bounds  on  target  loading  over  the  a;,  which 
define  7^  (E) ,  target  response  functions  having  specific  7g  (E)  were  chosen  as  test  cases. 
Resulting  optimized  bounds  and  7g(£)  are  compared  to  the  originally  defined  values. 

The  first  test  case  uses  a  linear  combination  of  instrument  response  functions  to  define 
T(E)  such  that  the  remainder  is  zero.  In  test  case  2,  T(E)  is  defined  as  the  same  linear 
combination  of  response  functions  but  T$  (E)  is  defined  as  a  positive  correction  outside 
the  instrument  response  envelope.  The  third  test  case  is  similar  to  the  preceding  case 
except  Tfi  (E)  is  defined  as  a  negative  correction;  the  same  T-  (E)  is  subtracted  from  the 

linear  combination  instead  of  being  added  to  it. 

Based  on  validation  results,  four  variations  were  performed  on  each  test  case.  The 
first  variation  defined  T(E)  as  a  positive  combination  of  closed  instrument  response 
functions.  The  second  variation  used  a  positive  combination  of  open  response  functions. 
The  third  variation  used  positive  combinations  of  both  types  of  detection  systems,  the  final 
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variation  was  not  restricted  to  linear  combinations  of  positive  terms  (i.e.,  the  a;  can  be 
negative).  The  selection  of  four  combinations  of  instrument  responses  requires  four 
Tfi(E)  definitions.  This  stems  from  a  non-negativity  constraint  on  the  target  response 
function.  The  four  linear  combinations  of  instrument  responses  are  listed  below. 


Variation  A: 
Variation  B: 
Variation  C: 


Variation  D: 


T(E)=R1(E)  +  2.5  R3(E) 

T(E)  =  R9(E)  +  Rn(E) 

T{E)  =  R\{E)  +  2  5R3(E)  +  R9(E )  +  Rn(E) 


T(E)  =  2 


\.2R\(E)  +  0.75  R2(E)  -  1.75  R7(E) 

-  0.15  R8(E)  +  0  5Rn(E)  +  R12(E) 


(41a) 

(41b) 

(41c) 

(4  Id) 


Figures  7-10  show  T(E)  and  Figures  11-14  show  the  defined  T^(E)  for  the  test  cases. 


Figure  7.  Test  Target  Response  Function,  Variation  A 
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Figure  8.  Test  Target  Response  Function,  Variation  B 


0  5  10  15  20  25 

Energy,  E/EO 


Figure  9.  Test  Target  Response  Function,  Variation  C 
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Figure  10.  Test  Target  Response  Function,  Variation  D 


Figure  11.  Defined  Test  T^{E) ,  Variation  A 
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Figure  14.  Defined  Test  (E) ,  Variation  D 


Examination  of  Figures  11-14  show  that  the  7-  (E)  values  are  relatively  small.  This  is  due 
to  the  non-negativity  constraint  on  T(E)  (a  concern  when  7^  (E)  is  negative)  and  the  fact 

that  the  T(E)  responses  shown  in  Figures  7-10  are  relatively  small  in  themselves. 

Results  of  Test  Case  1  (Tl).  Tl.A  defined  the  target  response  function  as  Eq 
(41a)  with  no  correction  term.  The  results  of  Tl.A  are  summarized  in  Table  5. 

Table  5  shows  numerically  calculated  total  target  loading  bounds  of  30%  for  a 
target  response  function  exactly  defined  as  a  positive  combination  of  instrument  responses. 
This  supports  the  analytic  argument  from  Chapter  II  that  through-folding  error  simplifies 
to  measurement  error  if  T(E)  is  an  exact  combination  of  positive  (or  negative)  terms. 
However,  neither  the  upper  nor  lower  defined  bounds  are  found  to  be  exactly  15%.  The 
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difference  between  expected  and  calculated  results  are  due  to  additional  numeric  error 


introduced  via  the  numeric  approximation  of  the  instrument  signal  when  a;  is  not  zero. 

Optimization  results  are  also  shown  in  Table  5.  Both  maximum  and  minimum 
values  were  found  via  through-fold  of  closed  detector  system  signals.  This  is  expected  as 
T(E)  is  defined  using  only  closed  system  response  functions.  The  optimization  routines 
were  unable  to  find  the  known  minimum  bound  (30%);  however,  the  total  bound  was 
reduced  from  100%  to  38.5%.  The  large  optimized  minimum  bound  (relative  to  maximum 
bound)  is  attributable  to  a  relatively  large  negative  value  of  as  on  error  propagation.  A 
complete  discussion  of  error  propagation  is  found  in  Bevington(1969:Ch  4). 


Table  5.  Test  Results  for  Case  Tl. A 


ai,  Naive 
Approach 

Defined 

ai 

Optimized 

Max 

Optimized 

Min 

1  al 

0 

1 

0.1284 

0.3427 

a2 

0 

0 

0.1778 

u 

a3 

0 

2.5 

2.5494 

1.3688 

a4 

0 

0 

0 

1.0799 

a5 

0 

0 

0.0221 

-0.9621 

a6 

0 

0 

0.0193 

0.7827 

a7 

0 

0 

0.0158 

0.5918 

a8 

0 

0 

-0.0026 

0.821 

a9 

0 

0 

0 

0 

alO 

0 

0 

0 

0 

all 

0 

0 

0 

0 

al2 

0 

0 

0 

0 

Zv( a),  %  + 

50 

13.5 

17.1 

n/a 

Zv{ a),  %  - 

50 

16.5 

n/a 

21.4 

Total  Bound,  % 

100 

30 

38.5 
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Tl.B  proved  to  be  a  trivial  variation  on  Tl.A.  The  linear  combination  of  open 
response  functions  from  Eq  (41b)  was  exactly  reproduced  by  through-fold.  The  resulting 
total  bound  was  found  to  be  30%.  Again,  a  slight  error  due  to  signal  computation  was 
experienced.  Results  from  TLB  are  contained  in  Appendix  C. 

Tl.C  also  proved  to  be  trivial,  despite  the  target  response  being  defined  as  the 
exact  linear  combination  of  both  open  and  closed  response  functions  given  by  Eq  (41c). 
All  linear  coefficients  were  of  the  same  sign  resulting  in  the  total  error  bound  being 
minimized  to  30%  using  defined  values  for  a,.  As  found  in  validation,  the  upper  and  lower 
optimized  bounds  did  not  result  from  through-fold  of  the  same  type  of  detector  systems. 
The  overall  optimized  error  bound  was  35.1%— a  slight  improvement  over  the  simpler 
model  from  TLA.  Results  from  Tl.C  are  contained  in  Appendix  C. 

T1.D  not  only  utilized  both  types  of  detector  systems  but  defined  T(E)  as 
combinations  of  positive  and  negative  terms  (i.e.,  positive  and  negative  aO  given  by  Eq 
(41d).  The  remainder  term,  7^  ( E ) ,  remains  defined  as  zero.  The  results  of  T1  .D  are 
summarized  in  Table  6.  The  influence  of  the  negative  coefficients  on  the  error  bound  is 
significant  for  the  defined  selection  of  ai.  The  total  defined  error  bound  increased  from 
30%  to  91%.  Because  T^(E)  is  zero,  the  increase  is  attributable  to  numeric  integration 
approximations  of  the  signals  and  to  propagation  of  errors  associated  with  these  signals 
due  to  the  defined  a;  (some  being  negative,  other  being  positive). 
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Table  6.  Test  Results  for  Case  Tl.D 


aj,  Naive 
Approach 

Defined 

ai 

Optimized 

Max 

Optimized 

Min 

al 

0 

4.2 

0.829 

0.762 

a2 

0 

0 

0.00002 

1.3613 

a3 

0 

1.5 

1.6527 

0.0002 

a4 

0 

0 

0 

1.4865 

a5 

0 

0 

3.8724 

-0.0002 

a6 

0 

0 

-1.7982 

1.7565 

al 

0 

-3.5 

4.7251 

2.3024 

a8 

0 

0 

-5.3757 

-1.7722 

a9 

0 

-0.15 

0 

0 

0 

0 

0 

1 

0 

0 

al2 

0 

2 

0 

0 

Zv( a),  %  + 

50 

48.6 

12.2 

n/a 

Zu(a),  %  - 

50 

42.8 

n/a 

22.6 

Total  Bound,  % 

100 

91.4 

34.8 

Upper  and  lower  optimized  error  bounds  were  both  determined  via  through-fold  of 
closed  instrument  signals.  Total  optimized  bounds  were  less  than  35%.  This  is  consistent 
with  values  found  in  previous  trials. 

The  importance  of  this  variation  is  not  only  that  through-fold  can  achieve  tighter 
bounds  than  a  defined  set  of  a\  but  that  the  optimal  (E)  is  not  always  the  minimum. 

Figure  11  plots  the  value  of  T^(E)  calculated  from  a;  for  optimized  upper  bound.  T^(E) 


was  originally  defined  as  zero.  Figure  15  shows  that  the  optimized  T^{E)  reflect  the 
general  characteristics  (peaks  and  valleys)  of  the  target  response  function  (Figure  10). 
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Results  of  Test  Case  2  (T2).  The  four  variations  of  T2  were  designed  to 
investigate  the  effect  on  optimization  when  the  defined  target  loading  function  includes  a 
positive  contribution  from  T#  (E) .  Since  the  results  of  T2.A,  T2.B,  and  T2.C  were  similar 

to  each  other,  their  results  will  only  be  briefly  discussed.  Results  from  T2.A,  T2.B,  and 
T2.C  are  contained  in  Appendix  C. 

The  first  three  variations  of  T2  produced  similar  results.  In  all  cases,  naive  error 
was  verified  to  be  ±50%,  the  defined  error  increased  from  ±15%  (analytically  and  in  Tl) 
to  ±24%,  and  the  optimized  error  bounds  remained  comparable  to  those  of  Tl.  Defined 
error  increase  is  expected  as  T(E)  is  not  an  exact  combination  of  instrument  responses. 

The  additional  error  is  associated  with  propagation  of  the  error  associated  with  the  non- 
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optimal  T5  (E)  propagating  into  the  exact  spectrum  (see  explanation  of  Eq  22b).  Unlike 


the  results  for  the  defined  error  bound,  total  optimized  bounds  remained  constant  with  the 
addition  of  a  remainder  term  to  T(E).  Since  the  optimized  bounds  from  T1  have  the 
effects  of  error  propagation  included,  large  additional  propagation  errors  are  not  expected 
in  the  optimized  results  of  T2.  Some  additional  error  is  expected  with  new  selections  over 
ag  however,  the  optimization  of  a;  minimizes  additional  error  propagation. 

Figure  16  compares  T2.C’s  optimal  T^(E)  for  the  optimized  upper  bound  to  the 
defined  Tg  (E) .  Similar  results  were  obtained  for  maximum  and  minimum  bounds  in 
T2.A,  T2.B,  and  T2.C.  Since  all  defined  a*  were  positive,  the  optimum  T(i  (E)  should  be 

no  larger  than  the  defined  Ta-(E) .  Numeric  integration  from  1 J E®  to  40Zsj°  of  the  curves 
in  Figure  16  yield  optimum  T$(E)  =  0.17,  and  defined  TS(E)  =0.34. 


Legend 
Def.  Ta 
Opt.  Ta 


Figure  16.  Case  T2.C,  Defined  and  Optimized  (Upper  Bound)  7 5  ( E ) 


Test  case  T2.D  used  combinations  of  both  types  of  detector  response  systems  and 


I 

allowed  ai  to  be  combinations  of  positive  and  negative  numbers.  The  remainder  term  was 


defined  positive  over  all  energies.  T2.D  results  are  presented  in  Table  7. 


Table  7.  Test  Results  for  Case  T2.D 


ai,  Naive 
Approach 

Defined 

at 

Optimized 

Max 

Optimized 

Min 

al 

0 

4.2 

0.829 

a2 

0 

0 

0.6938 

■m 

a3 

0 

1.5 

1.606 

0.454 

a4 

0 

0 

0.998 

2.0913 

a5 

0 

0 

4.2905 

0.0001 

a6 

0 

0 

1.2544 

2.7684 

a7 

0 

-3.5 

4.1012 

1.433 

a8 

0 

-0.15 

-4.614 

-0.6655 

a9 

0 

0 

0 

0 

alO 

0 

0 

0 

0 

all 

0 

1 

0 

0 

al2 

0 

2 

0 

0 

Zu( a),  %  + 

50 

42.1 

13.3 

n/a 

Zv{ a),  %  - 

50 

37.8 

n/a 

21.9 

Total  Bound,  % 

100 

79.9 

35.2 

As  seen  in  test  case  Tl.D,  upon  introduction  of  negative  coefficients  in  the 
definition  of  the  T(E),  the  defined  relative  error  bound  more  than  doubled  that  of  the  cases 
having  7g(£)  =  0.  Optimized  upper  error  bound  has  remained  at  35%.  Figure  17  plots 


the  defined  and  optimized  (upper  bound)  (E)  for  case  T2.D. 


Defined  and  Optimized  Remainder,  T2.D 


Legend 
Def.  Ta 
'  Opt.  Ta 

0  5  10  15  20  25  30 

_ Energy,  EO/E _ 

Figure  17.  Case  T2.D,  Defined  and  Optimized  (Upper  Bound)  T^(E) 

As  with  Tl.D,  the  optimized  7g(£)  has  no  similarities  to  the  defined  T$(E) .  However, 
the  shape  of  T2D’s  optimized  T^(E)  is  the  same  as  in  T1.D;  small  variations  are  expected 
(compared  to  T1D)  because  T(E)  has  changed  by  a  small  amount,  T^(E) . 

Results  of  Test  Case  3  (T3).  Test  case  three  consisted  of  the  same  linear 
combinations  of  instrument  responses  as  in  T1  and  T2  except  a  negative  Ts  (E)  was 

included.  For  consistency,  the  same  7- (E)  defined  for  cases  T2  were  used,  however, 

they  were  subtracted  from  Eqs  (41)  instead  of  being  added.  The  results  from  test  case  T3 
parallel  those  of  T2.  The  only  notable  differences  are  that  defined  and  optimized  bounds 
increased.  The  results  of  the  trivial  cases  (T3.A,  T3.B,  and  T3.C)  are  contained  in 
Appendix  C  and  are  briefly  discussed  below. 
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In  all  cases,  naive  error  was  verified  to  be  ±50%,  the  defined  error  increased  from 


±24%  in  T2  to  ±50%  in  T3,  and  the  optimized  error  bounds  increased  from  ±17%  in  T2  to 
±21%  in  T3.  The  increase  in  defined  error  bound  is  again  expected  due  to  the  negative 
contribution  of  T^(E) ,  in  effect,  mimicking  a  large  negative  a*. 

Figure  18  plots  the  defined  and  optimized  (upper  bound)  remainders  for  T3.C. 
Similar  results  were  obtained  for  upper  and  lower  bounds  in  T3.A  and  T3.B.  In  the 
optimizer’s  attempt  to  represent  the  complete  target  response  as  a  set  of  linear  terms,  the 
total  area  under  the  optimized  T^(E)  curve  in  Figure  18  was  driven  negative.  The  results 
from  the  overall  negative  T~  (E)  on  error  propagation  are  akin  to  selection  of  additional 
negative  ai  values.  Hence  the  error  bounds  associated  with  an  overall  negative  (E)  are 
expected  to  increase  over  other  cases. 


Figure  18.  Case  T3.C,  Defined  and  Optimized  (Upper  Bound)  T^(E) 
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The  results  from  T3.D  are  shown  in  Table  8.  The  overall  influence  of  T(E)  being 


explicitly  defined  with  a  negative  correction  term  is  an  increase  in  through-fold  bounds  to 
47%.  However,  the  bounds  are  still  a  considerable  improvement  over  the  naive  approach. 


Table  8.  Test  Results  for  Case  T3.D 


ai,  Naive 
Approach 

Defined 

ai 

Optimized 

Max 

Optimized 

Min 

al 

0 

4.2 

0.8284 

0 

a2 

0 

0 

0.5878 

0 

a3 

0 

1.5 

0.0871 

0 

a4 

0 

0 

0 

0 

a5 

0 

0 

2.6229 

0 

a6 

0 

0 

-1.4672 

0 

a7 

0 

-3.5 

4.679 

0 

a8 

0 

0 

-5.2097 

0 

a9 

0 

-0.3 

0 

0.1583 

alO 

0 

0 

0 

0.193 

all 

0 

1 

0 

0.6819 

al2 

0 

2 

0 

1.238 

Zu( a),  %  + 

50 

99.8 

19.1 

n/a 

Zv( a),  %  - 

50 

90.1 

n/a 

27.5 

Total  Bound,  % 

100 

189.9 

46.6 
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V.  Summary  and  Conclusions 

Summary 

Validation  of  Method.  The  purpose  of  the  validation  case  is  to  ensure  that  the 
through-fold  routine  is  operating  properly.  The  validation  case  was  set  up  with  two  goals: 
1)  verifying  through-folds  capabilities  to  simulate  the  case  of  having  no  experiment  and  2) 
verification  that  target  loading  bounds  can  be  decreased  from  the  no  experiment  case. 
Through-fold  proved  capable  of  predicting  pre-test  error  bounds  and  of  utilizing  test  data 
to  reduce  the  pre-test  prediction. 

Test  Cases.  The  overall  goal  of  the  test  cases  was  to  explore  the  feasibility  of 
through-fold.  A  series  of  test  cases  was  set  up  to  compare  the  through-fold  results  to  the 
expectations  for  the  method.  The  test  cases  were  set  up  to  incrementally  increase  the 
complexity  of  a  problem.  Test  case  one  assumed  target  loading  was  a  linear  combination 
of  instrument  responses.  Test  case  two,  increased  the  complexity  of  the  problem  by 
assuming  there  was  a  positive  remainder  term.  Finally,  test  case  three  allowed  the 
remainder  to  be  negative.  To  aid  in  evaluation,  each  test  case  consisted  of  four  separate 
combinations  of  detection  system  signals. 

All  basic  expectations  for  the  through-fold  were  met.  Results  from  test  case  one 
indicated  that  through-fold  could  predict  reasonable  bounds  on  an  unknown  function. 
Results  from  case  two  demonstrated  the  anticipated  behavior  concerning  error 
propagation.  The  results  from  the  more  complicated  problem  presented  in  test  case  three 
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showed  a  reduction  in  predicted  error  bounds  from  the  naive  case  and  from  the  error 
bound  based  on  the  defined  values  for  the  target  loading  response  functions. 

Conclusions 

The  results  from  the  validation  and  test  cases  support  the  following  conclusions: 

1.  Through-fold  consistently  reduced  the  error  bounds  associated  with  the  test 
cases.  However,  the  combination  of  energy  bin  resolution  and  small  spectrum 
and  target  loading  values  in  some  of  the  bins  may  have  limited  the  capabilities  of 
the  through-fold  routine. 

2.  Optimization  of  error  bounds  has  a  limiting  factor  dependent  on  propagation  of 
errors  and  the  actual  spectrum  and  response  functions  used. 

3.  The  optimization  does  not  always  occur  when  the  target  response  function  is 
defined  exactly. 

4.  Although  there  is  some  degree  of  user  judgment  involved  in  selection  of  the 
starting  values  for  the  optimization  routines,  the  resulting  target  loading  bounds 
did  not  appear  to  be  sensitive  to  the  selection. 

5.  The  simplified  target  loading  problem  was  adequately  handled  by  through-fold. 
The  results  of  test  cases  two  and  three  (where  through-fold  consistently 
outperformed  results  calculated  from  exact  target  response  functions)  indicate 
promise  for  the  technique. 
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VI.  Recommendation  s 


The  following  are  recommendations  for  continued  work  in  this  area.  First,  a  more 
robust  optimization  routine  must  be  found.  This  could  solve  a  few  of  the  problems 
experienced.  A  different  optimization  routine  may  speed  up  calculations,  should  be 
capable  of  finding  exact  solutions  for  problems  used  in  test  case  one,  and  should  be 
capable  of  handling  the  quadratic  problem  which  was  simplified  to  a  linear  one  for  this 
study.  It  is  hoped  that  insight  gained  into  the  through-fold  methodology  from  the  linear 
programming  problem  will  carry  over  to  quadratic  programming  problem.  Second, 
different  choices  of  input  data  should  be  used  to  verify  that  through-fold  results  are  similar 
in  more  realistic  test  cases.  Finally,  a  third  area  of  study  must  consider  how  well  through¬ 
fold  agrees  with  output  from  accepted  unfolding  routines. 
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Appendix  A:  Formal  Derivation  of  Through-Fold  Equations 


Definition  of  Terms 


5 


Ri(E)  = 


S(E)  = 
T(E)  = 

Yi 

Z 


tilde,  denotes  an  exact  quantity,  as  opposed  to  measured 

delta,  denotes  correction  term  determined  from  measurement  error, 
calibration  error,  or  a  conservative  estimate 

instrument  response  as  a  function  of  energy,  ith  instrument 

spectrum  as  a  function  of  energy 

target  response  as  a  function  of  energy 

instrument  signal,  ith  instrument 

target  loading 


NOTE:  The  absence  of  both  tilde  and  delta  indicates  the  quantity  is 

measured  or  calibrated. 


The  derivation  begins  with  a  definition  of  target  loading.  Target  loading  can  be 
viewed  as  the  signal  which  would  have  been  recorded  by  the  target  had  it  actually  been  an 
error  free  instrument.  Target  loading  can  be  represented  similarly  to  instrument  signal  by: 


Z  =  jS(E)T(E)dE  (Al) 

0 


None  of  the  values  in  Eq  (Al)  are  actually  known.  They  can,  however,  be  approximated 
by  a  measured  or  predicted  value  plus  a  correction  term.  If  the  predicted  target  response 
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is  further  approximated  as  a  linear  combination  of  the  instrument  responses  plus  a 
remainder,  to  keep  the  relationship  exact,  then  all  the  terms  on  the  right  hand  side  of  Eq 
(Al)  can  be  replaced  with  known  predicted,  or  boundable  values.  Defining  the  terms: 


S(E)  =  S(E)  +  SS(E) 

(A2.a) 

T(E)  =  T(E)  +  8T(E) 

(Alb) 

Ri(E)  =  Ri(E)  +  8Ri(E) 

(A2.c) 

Yi  =  Yi+6Yi 

(A2.d) 

Equations  (A2)  should  be  read  as  ‘the  exact  value  equals  the  measured  value  plus  a 
correction  term’.  The  four  correction  terms  in  Eqs  (A2)  are  based  on  a  conservative 
estimates  (the  first  two),  calibration  error,  and  measurement  error  respectively. 

Derivation 

Since  all  terms  on  the  right  hand  sides  (RHS)  of  Eqs  (A2)  are  considered  known, 
they  can  be  substituted  into  the  exact  target  loading  equation,  eliminating  the  unknowns. 
Substituting  Eq  (A2.b)  into  Eq  (Al)  yields: 

oo 

Z  =  |S(£)  [T(E)  +  ST(E)]  dE  (A3) 

0 

Now  let 

Z  =  Z  +  8Z  (A4) 
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where, 


(A4.a) 


Z=  f  S(E )  T(E)  dE 
0 


8Z=  js(E)  8T(E)  dE  (A4.b) 

0 


Assuming  T(E)  can  be  approximated  by  a  linear  combination  of  instrument 
response  functions  plus  a  correction  term  (remainder),  yields  the  following  approximation: 

T(£)  =  X«i  Ri(E)  +  Ta(£)  (A5) 

i 


where 


at  =  a  constant  for  the  ith  instrument 

Ta  (E)  =  a  remainder  based  on  the  choice  of  ai 


Substituting  this  last  equation  into  Eq  (A4.a), 


oo 


Z  =  \S{E) 
0 


'£aiRi(E)  +  Ta(E) 


dE 


(A6) 


or,  by  expanding  the  integral  and  reversing  the  order  of  integration  and  summation: 


Z  =  ^aijs(E)Ri(E)dE  +  \S(E)Td(E)dE  (A7) 

i  0  0 

The  exact  spectrum  in  Eq  (A7)  can  now  be  replaced  with  its  approximation  from  Eq 
(A2.a).  However  instead  of  replacing  S  (E)  in  both  integrals,  it  is  beneficial  to  replace 
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Ri(E)  in  the  first  integral  of  Eq  (A7)  with  its  value  from  Eq  (A2.c)  and  then  to  substitute 
the  approximation  of  the  exact  spectrum  into  the  second  part  of  Eq  (A7).  This  trick 
yields: 


oo  oo 

Z  =  YJai\S(E)Ri{E)dE  -  £aj5(E)5E2(E)dE 
i  0  i  0 

(A8) 

oo  oo 

+  J  S(E)  TS(E)  dE  +  J 85(E)  T5(E)  dE 
0  0 


The  above  trick  has  introduced  the  exact  signal,  the  first  integral  of  Eq  (A8),  into  our 
equation.  We  will  now  introduce  data  (measured  signals)  into  our  problem  by 
approximating  of  the  exact  signal  using  Eq  (A2.d). 


Z  =  ZaiYi  +  'ZaijS(E)bRi(E)dE 

i  i  i  0 

(A9) 

oo  oo 

+  J  S(E)Ts(E)dE  +  J5 S(E)Ta(E)dE 
0  0 


Substituting  the  approximation  of  the  exact  spectrum  for  S  ( E ) . 


z=  2"/1/  +  2>i5Ii 

i  i 

oo  oo 

-  ^fl/j5(E)5E/(E)dE  -  ^  aj  55(E)  8Ej(E)dE  (A10) 

i  0  /  0 

oo  oo 

+  J  5(E)  TS(E)  dE  +  J 55(E)  7g(E)  dE 
0  0 
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Now  that  all  the  exact  terms  have  been  eliminated  from  Z,  we  can  switch  our 
attention  to  the  8Z  term  from  Eq  (A4.b).  For  this  term  all  we  have  to  do  is  substitute  the 
approximation  of  the  exact  spectrum  and  simplify  it  to: 


oo  oo 

5 z  =  j  S(E)  8T(E)  dE  +  jbS(E)  bT(E)  dE  (All) 

0  0 


Finally,  we  can  combine  equations  (A10)  and  (All)  for  the  exact  target  loading: 

z=  2^  + 

i  i 

oo 

-  2>J5(£)  SRi(E)dE 
i  0 

OO 

+  \S(E)Ta(E)dE  + 

0 
oo 

+  J  S(E)  8T(E)  dE  + 


The  final  equation  for  target  loading  is  found  by  writing  the  above  equation 
explicitly  in  terms  of  the  constants,  ai,  by  replicating  T^(E)  with  its  value  from  Eq  (A5): 


-  '^ai\8S(E)dRi(E)dE 
i  0 

OO 

jbS(E)Ta(E)dE 


(M2) 


J  6S(E)  bT(E)dE 
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z=  +  5>,-  5Fi 


-  '^aijS(E)bRi(E)dE  -  ]TaJ  55(E)  5Rf(E)dE 


i  0 


I  0 


+  J  5(E)  T{E)~YjaiRi{E)  dE  +  Js 5(E)  T(E)-^aiRi(E)  dE 


(A  13) 


+  J  S(E)  §T(E)  dE  +  J  55(E)  bT(E)  dE 


Expanding  the  two  integrals  that  were  changed, 


z=  5>rf  +  J,ai&ri 


-  £a,JS,S(E)SR;(£)iffi 


i  0 


I  0 


+  J5(E)  7(E)  dE  -  J5(E)  ]Tfl.Rt(E)  dE 


0  j 


+  J  55(E)  7(E)  dE  -  J  55(E)  ]£«,•/?, -(E)  dE 
0  0  j 

OO  OO 

+  J  5(E)  57(E)  dE  +  J  55(E)  57(E)  dE 


(A  14) 


Reordering  the  terms,  to  include  reversing  the  order  of  integration  and  summation: 
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z=  Y.OiYi  +  X«i85i 

i  i 

oo  oo 

-  dE  -  5)  flj  J  8S(E)  E,(E)  dE 

i  0  i  0 

OO  oo 

-  Yaijs(E)5Ri(E)dE  -  Yaij&S(E)  hRi{E)dE  (A15) 

i  0  i  0 

OO  oo 

+  J  S(jB)  7(E)  dE  +  J  8S(£)  7(E)  dE 
0  0 

oo  oo 

+  jS(E)ST(E)dE  +  jSS(E)6T(E)dE 

0  0 

Now  that  we  have  an  expression  for  the  target  loading,  and  bounds  on  all  the  delta 
terms,  we  need  to  convert  the  analytic  expression  to  a  numeric  one. 


Conversion  to  a  Numeric  Solution 

There  are  two  forms  of  integrals  in  Eq  (A15)  which  need  conversion  to  numeric 
techniques,  those  with  the  instrument  response  function  involved  and  those  with  the  target 
response  function  involved.  Both  integrals  can  be  numerically  approximated  using  a 
summation  over  the  energy  range  of  concern.  Using  this  form  of  approximation,  basic 
instrument  signal  and  target  response  integrals  can  be  written  as: 


N 

yi  =  ^S(Ej)Ri(Ej)AEj  (A16.a) 

j= 1 

and 
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(A16.b) 


N 

z=y£S(Ej)T(Ej)AEj 

7=1 


where 

y>i  =  a  basic  approximation  to  signal  recorded  by  the  ith  instrument 

S(Ej)  =  the  predicted  spectrum  at  energy  step  j 

Ri(Ej)  =  the  response  function  of  the  ith  instrument  (as  measured) 

th 

A Ej  =  the  width  of  the  j  energy  step 

z  =  a  basic  approximation  to  target  signal 

T(Ej)  =  the  predicted  target  response  at  energy  step  j 


N  =  the  number  of  energy  groups 


Replacing  the  integrals  in  Eq  (A15)  with  summations  yields: 


Z=  5>,J;  +  2>,5y;- 


N  N 

-  '£ai'2dS(Ej)Ri(Ej)Mj  -  Jjai'£bS(Ej)  R^Ej)  AEj 

i  7=1  i  7=1 

N  N 

-  '£ai'£S(Ej)SRi(Ej)AEj  -  y£ai'2t&S{Ej)8Ri(Ej)  AEj 

i  7=1  i  7=1 

N  N 

+  ^S(Ej)T(Ej)AEj  +  ^6S(Ej)T(Ej)AEj 
7=1  7=1 

N  N 

+  ^SiEj)  bT(Ej)  AEj  +  £5S(E7)  ST(Ej)  AEj 

7=1  7=1 


(A  17) 
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To  make  the  above  calculation  easier,  it  is  beneficial  to  form  an  instrument  response 
matrix,  R,  a  spectrum  vector,  s,  and  a  target  response  vector,  t,  of  the  predicted  functions. 

(R)y  =Rl(E])AEJ 

(s  )j=S(Ej)  (A  18) 


(t)j=T(Ej)AEj 


Equations  (A  16 )  can  now  be  converted  to  vector  notation,  given  by 


N 

JjS(Ej)Ri(Ej)AEj  =  R  s 
7=1 


and 


(A  19) 


N 

^S(Ej)T(Ej)AEj=s-t 

7=1 


Utilizing  the  vector  notation  from  Eq  (A  18)  to  represent  the  summations: 


Z=  5>;Ji  + 

i  i 

-  ^  Oj  R-s  -  ^  ci[  R  •  5  s 

i  i 

-  ^  ai  5R  •  s  -  flj  8R  •  5s  (A20) 

i  i 

+  st  +  5s  t 
+  s-8t  +  5s-5t 
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Realizing  that  the  instrument  signals  can  be  represented  by  the  vector  (y);  =  [Yj]  and  that 
the  ai  can  be  represented  by  a  similar  vector  equation,  allows  the  previous  equation  to  be 
written  completely  in  terms  of  vectors  and  matrixes: 

Z  =  ay+  a-8y 

-  a  •  R  s  -  a  •  R  •  8s 

-  a  8R  s  -  a-8R-8s  (A21) 

+  s  t  +  8s • t 

+  s • St  +  8s -8t 


Making  use  of  the  distributive  property  of  dot  products  simplifies  the  exact  target  loading 
equation  becomes: 


Z  =  a  ■  (y  +  8y  -  R  s  -  R  •  8s  -  8R  •  s  -  8R  ■  8s) 
+  s  t  +  8s  t  +  s-St  +  8s -8t 


(A22) 


Equation  (A22)  is  an  exact  expression  for  target  loading  because  the  exact  correction 
vectors,  8s,  8y,  8R,  and  8t  have  been  carried  throughout  the  derivation.  Unfortunately, 

Eq  (A22)  can  not  be  computed  because  the  exact  delta  values  are  unknown. 

Bounds  on  Exact  Target  Loading.  The  value  of  Z  given  by  Eq  (A22)  cannot  be 
calculated  because  the  delta  terms  are  not  known  exactly,  only  within  the  bounds  on  the 
correction  terms  discussed  previously.  The  exact  target  loading  can  be  bound  by 

ZL  <Z<ZU  (A23) 
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where  Zu  and  ZL  are  the  upper  and  lower  bounds  respectively  on  target  loading. 

The  upper  (lower)  target  loading  bounds  defined  above  can  be  computed  for  any 
fixed  choice  of  a  by  choosing  the  values  for  5s,  8y,  8R,  and  8y  —  all  being  allowed  to  vary 
(within  the  error  bounds)  at  each  energy  bin  --  to  maximize  (minimize)  the  target  loading 
expression,  which  is  implicitly  a  function  of  a.  These  bounds  are  computed  by  solving  the 
following  linear  programs  for  the  upper  and  lower  target  loading  bounds  respectively: 


ZU(  a) 


=  Max 
8y,8s,8R 


a  •  (y  +  8y 
+  s  t 
,St 


-  R  s  -  R  5s 
+  8s  t  +  s-8t 


-  8R  s  - 
+  8s -8t 


8R  •  5s) 


(A24) 


ZL(  a) 


Min 


a(y 


+ 

+ 


8y,Ss,8R,8t 


8y  -  Rs 
st  +  8st 


R  •  8s  -  8R  •  s  - 
+  s-8t  +  5s -5t 


8R  •  8s) 


In  Eq  (A24)  8s,  8y,  8R,  and  8y  are  the  variables  subject  to  the  following  constraints: 

^min  —  ^y  —  5y  max 
^smin  —  Ss  <  5smax 

(A25) 

^^min  —  —  ^^max 

^min  ^  8t  <  5tmax 


R,  s,  t,  and  y  are  input  data;  and  a  is  the  variable  for  which  the  bound  is  to  be  determined. 
The  goal  functions  of  both  bounds  in  Eq  (A24)  are  the  same  and  are  given  by: 
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/(Sy,8s,8R,St)  =  a(y  +  Sy  -  Rs  -  R8s  -  SRs  -  SR -8s) 

(A26) 

+  s  t  +  8s  t  +  s  •  8t  +  Ss-8t 

The  upper  bound  is  found  by  maximum  the  goal  function  over  the  8s,  8y,  8R,  and  8y. 
Similarly,  the  lower  bound  is  found  by  minimizing  the  goal  function  over  the  8s,  8y,  8R, 
and  8y.  The  resulting  values  for  Zu  ( a)  and  ZL( a)  from  Eq  (A24)  represent  the  largest 
possible  bounds  (maximum  and  minimum  respectively)  on  target  loading  for  fixed  set  of  a. 
However,  it  should  be  possible  to  tighten  these  bounds  by  selecting  appropriate  sets  of  a. 
The  process  of  selecting  the  a’s  resulting  in  the  tightest  possible  bounds  is  accomplished 
by  optimizing  the  upper  and  lower  bound  over  all  possible  choices  of  a.  The  result  is  a 
classic  minimum-maximum  (mini-max)  problem  given  by: 


4nin  =  Mm[zy(a)] 
a 

(A27) 

Zmax=Max[zL(a) 

a 


Z^n  is  the  minimum  upper  bound  (least  upper  bound)  for  target  loading  and  is  found  by 
minimizing  Zu  (a)  over  a.  Z^ax  is  the  maximum  lower  bound  for  target  loading,  found 
by  maximizing  ZL(a)  over  a.  The  outer  loop  of  the  mini-max  problem  produces  an 
unconstrained  (the  a  vector  can  take  on  any  values)  quadratic  problem. 
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Standard  Form.  The  goal  function  in  Eq  (A26)  can  not  directly  be  used  in  a 
quadratic  programming  problem  because  it  is  not  in  standard  form.  The  standard  form  for 
linear  programming  problems  generally  require  non-negative  variables.  The  bounded  error 
terms  are  therefore  defined  as  the  difference  of  two  positive  numbers 

Sy  =  5y+  -  8y- 
8s  =  8s+  -  8s 

(A28) 

SR  =  8R+  -  8R 
St  =  8t+  -  8t 


where  all  the  terms  on  the  RHS  of  Eq  (A28)  are  greater  than  or  equal  to  zero. 

The  inner  loops  of  the  mini-max  problems  will  force  one  of  the  terms  on  the  RHS 
of  Eq  (A28)  to  zero  and  the  other  to  its  maximum  magnitude  at  each  energy  bin.  The  first 
term  will  be  zero  if  the  error  is  negative  while  the  second  term  will  be  zero  if  the  error  is 
positive.  Inserting  the  non-negativity  constraints  into  the  goal  function,  Eq  (A26): 


y 


/(8y,8s,8R,8t)  =  a-^ 


+  (8y+  -  8y  )  -  R  s 
R  (8s+  -  8s  )  -  (8R+  -  8R)s  • 
(8R+  -  SR-)(8s+-  8s-)] 


(A29) 


+  st  +  (8s +  -  8s  )  t  +  s  (8t+-  8t  ) 
+  j(8s+ -  8s  )  (8t+-  8t  )] 


or 
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y  +  Sy+  -  8y  -  R  •  s 

/(Sy,8s,8R,8t)  =  a  —  R  8s+  +  R-  8s  -  8R+s  +  8R~s 

-  [8R+  -Ss+  -  SR+-8s~  -  8R  -8s+  +  SR~Ss~l 
l  (A30) 

+  s  t  +  8s+t  -  8s  t  +  s  8t+  -  s-8t 
+  8s+St+  -  8s+-8t  -  8s  -8t+  +  8s  • St  J 

This  final  goal  function  remains  quadratic  in  the  delta  terms  (as  is  the  analytic  solution) 
and  is  subject  to  the  following  constraints: 

0<8y+  <8ymax 
0<8y“  <|-5ymin| 

0<8s+  <8smax 
0<8s  ^j-5smin| 

(A31) 

0<8R+  <8Rmax 
0<8R_  <|-8Rmin| 

0<8t+<Stmax 

0<8r<|-8tmin| 

The  quadratic  goal  function  Eq  (A30)  can  be  simplified  to  the  linear  goal  function  from 
Chapter  2  by  assuming  8R  and  8t  are  negligible  compared  to  8y  and  8s  and  setting  them 
equal  to  zero. 
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Appendix  B:  Final  Mathematica  Code 


A  Mathematica  notebook  for  through-fold 
of  bounds  on  exact  target  loading 


SECTION  L 

Data  Input  and  Definitions  of  Functions 

Enter  the  following  data  into  the  Data  Input  block  then  evaluate  the  cell: 


Energy  range  required:  Erange 

number  of  bins  covering  eRange:  numBins 

number  of  closed  instruments:  nClosed 

number  of  open  instruments:  nOpen 

fluorescer  k-edge  of  1  st  detector:  eOinitial 

bin  spacing  (geo  or  lin):  spacing 


eRange=64; 

nClosed=8; 

nOpen=4; 

numBins=30; 

spacing=geo; 

Spectrum  Definition: 

Choose  a  spectrum  or  create  your  own.  Enter  necessary  information 
and  evaluate  cells.  Evaluate  plot  cell  to  view  your  selection. 

Note:  only  one  spectrum  can  be  selected  at  a  time.  Evaluate  the  ClearAll[s  ] 
cell  prior  to  defining  a  second  spectrum. 

Clear All[s] 

Planckian  (1  or  2  temp) 

Spectrum  Data: 
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kTl=l; 

scalel=l; 

kT2=l; 

scale2=0; 

Spectrum: 

s[0]=0; 

s[0.0]=0; 

s[e J :=N[scalel ( 1 5/(kTl  Pi)A4)(eA3/(Exp[e/kTl]-l))  + 
scale2(  1 5/(kT2  Pi)A4)(eA3/(Exp[e/kT2]- 1 ))] ; 

Plot[s[e] ,  { e,0,eRange } , 

PlotLabel->"Planckian  Spectra", 

Frame->True, 

FrameLabel->{ "Energy,  E/EO",  "Amplitude"}, 
PlotRange->All]; 


Target  Response  Definition: 

The  same  guidance  applies  as  for  selection  of  the  spectrum. 
Clear  All  [t] 


Broad  Gaussian  ( change  shape  via  mean  and  standard  deviation) 
Gaussian  Data: 


mean=0; 

sd=15; 

sqrt2pi=N[Sqrt[2Pi]] ; 

Gausian  Response: 

t[e_]:=(l/(sd  sqrt2pi))  Exp[(-l/2)((e-mean)/sd)A2] 

PlotLabel->"Gausian  Response  Function", 
Frame->True, 

FrameLabel->{ "Energy,  E/EO",  "Amplitude"}, 
PlotRange->All]; 
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Detector  Systems  Definition: 

Evaluate  the  cells  associated  with  this  section.  They  are  based  on  previous  input. 


K-Edges 

numInst=nClosed  +  nOpen; 
sqrt2=N[Sqrt[2]]; 

e0closed=Table[e0initial(sqrt2An),  { n,0,nClosed- 1 }] ; 
elclosed=Table[2  e0initial(sqrt2An),  {n,0,nClosed-l}]; 
eOopen=Table[eOinitial  (2An),  { n, 0, nOpen- 1}]; 
e0all=Join[e0closed,  eOopen]; 


Detector  Response  Functions 

rclosed[eOclosed_,elclosed_]  := 

Which  [e<e0closed,0, 

e()closed<=e<e  1  closed, 
(l/e)(l-Exp[-2  (e0closed/e)A3]) 
(Exp  [(- 1/4)  (e  1  closed/e)  A3] ), 
elclosed<=e, 

(l/e)(l-Exp[-2  (e0closed/e)A3]) 
(Exp[(-3/2)(elclosed/e)A3]) 

] 

ropen[eOopen_]  :=If[e<e0open,0, 

(l/e)(l-Exp[-3  (e0open/e)A3]) 

]; 


Instrument  Responses 

Do[r[i,e_]=rclosed[eOclosed[[i]],elclosed[[i]]], 

{i,nClosed} 

]; 

Do[r  [i+nClosed,e_]  =ropen[eOopen[  [i]]] ,  { i, nOpen }  ] ; 
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Energy  Bins  Definition: 

Evaluate  the  cells  associated  with  this  section.  They  are  base  on  previous  input. 


Bin  Spacing: 

flag=spacing; 

deltaEgeo:=Table[N[eOinitial(eRangeA(i/numBins)  - 
eRangeA((i-  l)/numBins)) 

], 

{i,l,numBins} 

]; 


deltaElin:=N[Table[(eRange- 1)  eOinitial/numBins, 
{numBins}]]; 

dE=If[flag==geo,  deltaEgeo,  flag==lin,  deltaElin]; 


Bin  Boundaries 

energyBoundary  [0]  =eOinitial; 

eBound=Prepend[Table[energyBoundary[i]  =  dE[[i]]  + 
energyBoundary[i- 1  ] ,  { i, numBins }  ]  ,eOinitial]; 

eLower=Drop[eBound,- 1]; 
eUpper=Rest[eBound] ; 


Bin  Center  Energies  (for  evaluation  of  functions) 

eGeoCen=Sqrt[eUpper  eLower]; 
eLinCen=(eUpper  +  eLower)/2; 

eCen=If[flag==geo,  eGeoCen,  flag==lin,  eLinCen]; 
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SECTION  II. 


Evaluation  of  Vectors  and  Matrices 

Evaluate  the  cells  associated  with  this  section.  They  are  based  on  previous  input. 

Evaluate  the  Spectrum  Vector 

sVec=Table[s[eCen[[i]]] ,  { i,numBins }  ]; 

Evaluate  the  Target  Response  Vector 

tCen=Table[t[eCen[[i]]],{i,numBins}]; 
tVec=tCen  dE 


Evaluate  the  Instrument  Response  Matrix 

resmatrix=Table[(r[i,#]&  /@  eCen)  dE, 
{i,numlnst} 

]; 


Evaluate  the  Instrument  Signal  Vector 
yCen=resmatrix.sVec 


Add  Random  Noise  to  Signal  (Note:  Ym=Y+noise[  ]) 

Noise  should  NOT  be  added  if  using  assuming  no  measurement  error. 

If  this  is  the  case,  define  yM=yCen  +  0  prior  to  evaluation  of  the  yM  cell. 

«Statistics'  ContinuousDistributions' ; 

SeedRandom[l]; 

noise[frac_]:=Table[Random[NormalDistribution[0,frac]], 

{i,numlnst}]; 

measErr=noise[ .  15]; 

yM=yCen  +  (measErr*yCen) 
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Define  All  Delta  Terms 


The  delta  terms  are  set  for  up  to  60  bins,  any  spacing,  and  up  to  17  detectors. 

For  addiional  bins  and  instruments,  the  definitions  must  be  extended. 

The  form  is  self  evident. 

If  assuming  no  measurement  error,  dYmax  and  dYmin  should  be  zero  vectors. 

dSpositive=  { dSp  1 ,dSp2,dSp3,dSp4,dSp5,dSp6,dSp7,dSp8, 
dSp9,dSplO,dSpl  l,dSpl2,dSpl3,dSpl4,dSpl5, 
dSp  16,dSp  17,dSp  1 8,dSp  19,dSp20,dSp2  j  ,dSp22, 
dSp23,dSp24,dSp25,dSp26,dSp27,dSp28,dSp29, 
dSp30,dSp3 1  ,dSp32,dSp33,dSp34,dSp35,dSp36, 

dSp37,dSp38,dSp39,dSp40,dSp41,dSp42,dSp43, 

dSp44,dSp45,dSp46,dSp47,dSp48,dSp49,dSp50, 
dSp5 1  ,dSp52,dSp53,dSp54,dSp55,dSp56,dSp57, 
dSp58,dSp59,dSp60}; 

dSnegative={dSnl,dSn2,dSn3,dSn4,dSn5,dSn6,dSn7,dSn8, 
dSn9,dSn  1 0,dSn  1 1  ,dSn  1 2,dSn  1 3  ,dSn  1 4,dSn  1 5 , 
dSnl6,dSnl7,dSnl8,dSnl9,dSn20,dSn21,dSn22, 
dSn23,dSn24,dSn25,dSn26,dSn27,dSn28,dSn29, 
dSn30,dSn3 1  ,dSn32,dSn33,dSn34,dSn35,dSn36, 

dSn37,dSn38,dSn39,dSn40,dSn41,dSn42,dSn43, 

dSn44,dSn45,dSn46,dSn47,dSn48,dSn49,dSn50, 
dSn5 1  ,dSn52,dSn53,dSn54,dSn55,dSn56,dSn57, 
dSn5  8  ,dS  n59,dSn60 } ; 

dYpositive=  { dYp  1  ,dYp2,dYp3  ,dYp4,dYp5,dYp6,dYp7,dYp8 , 
dYp9,dYp  10,dYp  1 1  ,dYp  12, dYp  1 3, dYp  14, dYp  15, 
dYpl6,dYpl7}; 

dY  negative^  { dYn  1  ,dY  n2,dYn3  ,d  Y  n4  ,dYn5  ,dY  n6,dY  n7,dY  n8 , 
dYn9,dYnlO,dYnll,dYnl2,dYnl3,dYnl4,dYnl5, 
dYnl6,dYnl7}; 

dSmax=Table[0.5,  {60}]; 
dSmin=Table[0.5,  {60}]; 
dYmax=Table[0. 15,  {17}]; 
dY min=Table[0. 1 5,  {17}]; 
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Write  Vectors  to  Data  Files 


sVec»spectrum; 

tVec»target; 

resmatrix»rMatrix; 

yM»instSig; 

dY  positive»dY  pos; 
dY  negative»d  Y  neg; 
dYmax»dYplus; 
dYmin»dY  minus; 

dSpositive»dSpos; 
dS  negative»dS  neg; 
dSmax»dSplus; 
dSmin»dsminus; 


***  Additional  Information  {does  not  need  to  be  evaluated)  *** 


Exact  Instrument  Signal 

Do[yl[i]=NIntegrate[s[e]*r[i,e], 

{ e,eOall[[i]  ]  ,e  1  closed[  [i]],  64}], 
{i,nClosed}]; 

Do[yl[i]=NIntegrate[s[e]*r[i,e], 

{ e,eOall[[i]]  ,64 }  ] ,  { i,nClosed+ 1  ,numlnst }  ] ; 

yExact=Table[y  1  [i] ,  { i,numlnst } 


Exact  Target  Loading 

zExact=NIntegrate[s[e]  t[e],  {e,  0,  eRange}] 


Pseudo-Exact  Target  Loading 


pseudoZexact=s.t 


SECTION  in. 


Definition  and  Optimizationof  Exact  Target  Loading  Bounds 


-  Produces  dummy  files  from  datafiles  of  the  required  size. 

-  Defines  target  loading,  maximum,  and  minimum  functions. 

-  Provides  input  statements  for  calculation  of: 

-  Reads  full  information  from  datafiles. 

-  Results  for  no  experiment  (a's=0). 

—  FindMinimumjMAXX] 

—  FindMinimum[MINN] 

Notes:  1 )  File  written  to  and  read  from  the  Mathematica  default  directory. 

2)  Numbers  in  filenames  indicate  size  of  file  (inst  x  bins). 

3)  This  file  must  be  run  to  create  the  target  loading  equation  and  the 
expressions  for  use  in  the  maxx  and  minn  commands. 

4)  The  sections  starting  with  "How  Many  Instruments  and  Bins  in 
this  Test  ??"  and  ending  just  before  "Sample  Input  Statements 
for  Results"  must  be  evaluated  for  each  combination  of 
instruments  that  results  are  desired  for. 

Remember  to  change  filenames  appropriately. 


How  Many  Instruments  and  Bins  to  Through-Fold  ?? 

If  only  desiring  the  open  instruments,  enter  the  TOTAL  number  of  instruments 
(closed  +  open)  up  to  the  instrument  number  you  are  interested  in. 

inst=8; 

bins=30; 


Read  Information  From  Disk  Into  Dummy  Variables 

spec=«spectrum; 

targ=«target; 

resMat=«rMatrix; 

signal=«instSig; 
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dYpositive=«dYpos; 
dY  negative=«dY  neg; 
dY  max=«dY  plus; 
dY  min=«dY  minus; 

dSpositive=«dSpos; 
dS  negative=«dS  neg; 
dSmax=«dSplus; 
dSmin=«dsminus; 


Read  Appropriate  Parts  of  Dummy  Variables  Into  Vectors: 

If  only  want  open  instruments  and  the  vectors  include  both,  change  the 
instrument  counters  from  "{ i,inst }"  to  ”{i,nclosed+l  ,inst}"  for  R,  y, 
dYpos,  dY  neg,  dYplus,  and  dY minus  .  It  may  be  easiest  to  enter 
the  actual  number  instead  of"nclosed+l". 

s=Table[spec[[i]|,  {i,bins}]; 
t=T able [ tar g [ [i] ] ,  {i,bins}]; 

R=Table[resMat[[i,j]] ,  { i,inst } ,  { j ,bins } ]; 
y=Table[signal[[i]] ,  {i,inst}]; 

dSpos=Table[dSpositive[[i]] ,  { i,bins }  ] ; 
dSneg=Table[dSnegative[[i]],  {i,bins}]; 
dSplus=Table[dSmax[[i]],  { i,bins }]; 
dSminus=Table[dSmin[[i]],  { i,bins }]; 

dYpos=Table[dYpositive[[i]],  {i,inst}]; 
dY neg=T able [d Y negative [  [i]  ] ,  { i ,inst }  ] ; 
dYplus=Table[dYmax[[i]],  {i,inst}]; 
dY minus=T able[dYmin  [[i]  ] ,  { i,inst }  ] ; 


Checking  vectors  to  ensure  the  right  dimensions. 

Dimensions[y] 

Dimensions[R] 

Dimensions[s] 

Dimensionsft] 
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Definition  of  Functions 


Target  Loading  (z): 

terml=R.s; 
term2=R.(dSpos  s); 
term3=R.(dSneg  s); 
term4=(dSpos  s).t; 
tenn5=(dSneg  s).t; 

z[a_List,dYpos_List,dYneg_List,dSpos_List, 

dSneg_List]= 

a.(y  +  dYpos  y  -  dYneg  y  -  terml  - 
term2  +  term3)  + 
s.t  +  term4  -  term5; 

Writing  target  loading  funtion  (in  Mathematica  notation)  to  a  file 
and  verifying  it  is  there. 

Note:  Only  evaluate  the  second  command  if  you  WANT  to  see  a  huge  eqn. 
or  if  you  are  not  sure  the  file  was  written  to  disk  and  want  to  check. 

NOMENCLATURE:  File  names  indicate  target  loading  (z),  validation  (val), 

and  eight  instruments  in  the  equation  (8). 


%»z_val8; 

Hz_val8 


Maximum  Target  Loading  (maxx): 

This  is  a  DIFFERENT  correction  than  for  z. 

If  only  want  open  instruments  and  the  vectors  include  both,  change  the  instrument 
counter  from  "{ i,inst }"  to  "{ i,inst  -  nclosed}"  for  dYpos,  dYneg,  dYplus,  and  dYminus. 
It  may  be  easiest  to  enter  the  actual  number  instead  of  "inst  -  nclosed”. 

Note:  There  are  two  counters  for  dYpos  and  dYneg  each! 
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maxx[a_List]= 


ConstrainedMax[z[a, 

Evaluate[T  able[dY  pos  [[i]  ] ,  { i,inst }  ] , 
Table[dYneg[[i]],{i,inst}], 
Table[dSpos[[i]] ,  { i,bins }  ] , 
Table[dSneg[[i]],{i,bins}] 

] 

], 


Evaluate[Flatten[  { 

Table[dYpos[[i]]<=dYplus[[i]],  {i,mst}], 
Table[dYneg[[i]]<=dYminus[[i]],{i,inst}], 
Table[dSpos[[i]]<=dSplus[[i]],  {i,bins}], 
Table[dSneg[[i]]  <=dSminus[[i]] ,  { i,bins }  ] 
}] 

], 


Evaluate  [Flatten  [{ Table  [dYpos  [  [i]  ] ,  { i  ,inst }  ] , 
Table[dYneg[[i]],  {i,inst}], 
Table[dSpos[[i]],  {i,bins}], 
Table[dSneg[[i]],  {i,bins}] 

}] 


]; 


%»maval8; 


Minimum  Target  Loading  (minn): 

This  is  the  SAME  correction  as  for  minn. 

If  only  want  open  instruments  and  the  vectors  include  both,  change 
the  instrument  counter  from  "{ i,inst }"  to  "{i,inst  -  nclosed}"  for  dYpos, 
dYneg,  dYplus,  and  dYminus.  It  may  be  easiest  to  enter  the  actual 
number  instead  of  "inst  -  nclosed". 

Note:  There  are  two  counters  for  dYpos  and  dYneg  each! 

minn[a_List]= 
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ConstrainedMin[z[a, 
Evaluate[Table[dYpos[[i]],{i,inst}], 
Table[dY  neg[[i]] ,  { i,inst }  ], 
Table[dSpos[[i]] ,  { i,bins }  ], 
Table[dSneg[[i]],{i,bins}] 

] 

], 


Evaluate[Flatten[{ 

Table[dYpos[[i]]<=dYplus[[i]],  {i,inst}], 
Table[dYneg[[i]]<=dYminus[[i]],{i,inst}], 
Table[dSpos[[i]]<=dSplus[[i]],  {i,bins}], 
Table[dSneg[[i]]<=dSminus[[i]],{i,bins}] 
}] 

], 


Evaluate[Flatten[{Table[dYpos[[i]],  {i,inst}], 
Table[dYneg[[i]],  {i,inst}], 
Table[dSpos[[i]],  {i,bins}], 
Table[dSneg[[i]],  {i,bins}] 

}] 


]; 


%»mival8; 


Sample  Input  Statements  for  Results 

The  form  of  these  input  statements  will  be  consistent;  however,  there 
are  obvious  changes  based  on  the  number  of  instruments  used  and  the 
starting  points  desired. 


Read  datafiles  from  disk 
Target  Loading: 

z[a_List,dYpos_List,dYneg_List,dSpos_List, 
dS  neg_List] =«z_val8 ; 

Constrained  Maximum: 

maxx  [aJList] =«maval8 ; 
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Constrained  Minimum: 


minn[a_List] =«mival8 ; 

No  Experiment,  ai's  =  0. 

MAXX=maxx[  { 0,0, 0,0, 0,0, 0,0 }  ] 

MINN=minn[  { 0,0, 0,0, 0,0, 0,0 }  ] 

Optimizing  MAXX  and  MINN  over  the  Choice  of  ai's 

Minimizing  the  maximum. 

FindMinimum[First[maxx[  { al  ,a2,a3,a4,a5,a6,a7, 

a8 }]], 

{al,-l,l  },{a2,l,-l  },{a3,-l,l  },{a4,l,-l }, 

{a5,-l,l  },{a6, 1,-1  },{a7,-l,l  },{a8,l,-l }] 

Maximizing  the  minimum. 

Since  we  are  actually  minimizing  the  negative  of  the  minimum,  the  result  will 
be  the  NEGATIVE  of  the  true  answer. 

FindMinimum[-First[minn[{al,a2,a3,a4,a5,a6,a7,a8}]], 

{al,-l,l  },{a2,l,-l  },{a3,-l,l  },{a4,l,-l }, 

{a5,-l,l  },{a6, 1,-1  },{a7,-l,l  },{a8,l,-l }] 
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Appendix  C:  Additional  Test  Results 


This  appendix  is  a  continuation  of  the  results  of  Chapter  IV.  The  inputs  are  the 
same  as  defined  in  Chapter  IV. 


Table  9.  Test  Results  for  Case  Tl.B 


ai,  Naive 
Approach 

Defined 

at 

Optimized 

Max 

Optimized 

Min 

al 

0 

0 

0 

0 

a2 

0 

0 

0 

0 

a3 

0 

0 

0 

0 

a4 

0 

0 

0 

0 

a5 

0 

0 

0 

0 

a6 

0 

0 

0 

0 

a7 

0 

0 

0 

0 

a8 

0 

0 

0 

0 

a9 

0 

1 

1 

1 

alO 

0 

0 

0 

o  1 

all 

0 

1 

1 

1 

al2 

0 

0 

0 

0 

Zu( a),  %  + 

50 

18.6 

18.6 

n/a 

Zu( a),  %  - 

50 

12.3 

n/a 

12.3 

Total  Bound,  % 

100 

30.9 

30.9 
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Table  10.  Test  Results  for  Case  Tl.C 


&i.  Naive 

Defined 

Optimized 

Optimized 

Approach 

at 

Max 

Min 

al 

0 

1 

2.4452 

0 

a2 

0 

0 

1 

0 

a3 

0 

2.5 

2.0121 

0 

a4 

0 

0 

0.1745 

0 

a5 

0 

0 

3.6246 

0 

a6 

0 

0 

1.5946 

0  ! 

a7 

0 

0 

1.3151 

0 

a8 

0 

0 

-1.9238 

0 

a9 

0 

1 

0 

1.3333 

alO 

0 

0 

0 

0.923 

all 

0 

1 

0 

0.9613 

al2 

0 

0 

0 

0.035 

Zu(a),  %  + 

50 

15.8 

18 

n/a 

Zu( a),  %  - 

50 

14.4 

n/a 

17.1 

Total  Bound,  % 

100 

30.2 

35.1 
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Table  11.  Test  Results  for  Case  T2.A 


ai,  Naive 
Approach 

Defined 

ai 

Optimized 

Max 

Optimized 

Min 

al 

0 

1 

0.769 

0.4971 

a2 

0 

0 

0.1211 

0.9648 

a3 

0 

2.5 

3.3036 

2.5402 

a4 

0 

0 

0 

0.9931  ! 

a5 

0 

0 

-0.0001 

-1.1683 

a6 

0 

0 

-0.1288 

0.826 

a7 

0 

0 

0.0702 

-0.6397 

a8 

0 

0 

0.0897 

0.8914 

a9 

0 

0 

0 

0 

alO 

0 

0 

0 

0 

all 

0 

0 

0 

0 

al2 

0 

0 

0 

0 

Zu(a),  %  + 

50 

23.6 

18.5 

n/a 

Zu( a),  %  - 

50 

26 

n/a 

17.6 

Total  Bound,  % 

100 

49.6 

36.1 
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Table  12.  Test  Results  for  Case  T2.B 


ai,  Naive 
Approach 

Defined 

ai 

Optimized 

Max 

Optimized 

Min 

al 

0 

0 

2.0801 

0 

a2 

0 

0 

0.9988 

0 

a3 

0 

0 

0.6178 

0 

a4 

0 

0 

0.1997 

0 

a5 

0 

0 

5.2567 

0 

a6 

0 

0 

2.0032 

0 

a7 

0 

0 

1.2941 

0 

a8 

0 

0 

-2.2269 

0 

a9 

0 

1 

0 

1.1834 

alO 

0 

0 

0 

-0.03 

all 

0 

1 

0 

1.4374 

al2 

0 

0 

0 

-0.1321 

Zv (a),  %  + 

50 

26.7 

21.6 

n/a 

Zu(a),%~ 

50 

22 

n/a 

18 

Total  Bound,  % 

100 

48.7 

39.6 
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Table  13.  Test  Results  for  Case  T2.C 


ai,  Naive 
Approach 

Defined 

ai 

Optimized 

Max 

Optimized 

Min 

al 

0 

1 

2.4452 

0 

a2 

0 

0 

1 

0 

a3 

0 

2.5 

2.8914 

0 

a4 

0 

0 

0.612364 

0 

a5 

0 

0 

4.6113 

0 

a6 

0 

0 

-0.475 

0 

a7 

0 

0 

0.001 

0 

a8 

0 

0 

-0.6234 

0 

a9 

0 

1 

0 

1.3333 

alO 

0 

0 

0 

1.2265 

all 

0 

1 

0 

0.3256 

al2 

0 

0 

0 

-0.0881 

Zu (a),  %  + 

50 

22.4 

15.5 

n/a 

Zu(  a),%- 

50 

21.2 

n/a 

18.8 

Total  Bound,  % 

100 

43.6 

34.3 
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Table  14.  Test  Results  for  Case  T3.A 


a;,  Naive 
Approach 

IH 

Optimized 

Max 

Optimized 

Min 

al 

0 

1 

-0.3123 

0 

a2 

0 

0 

0.4192 

0 

a3 

0 

2.5 

1.4946 

0 

a4 

0 

0 

-0.4488 

0 

a5 

0 

0 

0.7568 

0 

a6 

0 

0 

0.459 

0 

a7 

0 

0 

0.3681 

0 

a8 

0 

0 

-0.5048 

0 

a9 

0 

0 

0 

0.0993 

alO 

0 

0 

0 

0.4922 

all 

0 

0 

0 

0.0711 

al2 

0 

0 

0 

0.0218 

Zv( a),  %  + 

50 

55.3 

25.6 

n/a 

Z^(a),  %  - 

50 

60.9 

n/a 

19.5 

Total  Bound,  % 

100 

116.2 

45.1 
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Table  15.  Test  Results  for  Case  T3.B 


ai,  Naive 
Approach 

Defined 

a; 

Optimized 

Max 

Optimized 

Min 

al 

0 

0 

0.4492 

0 

a2 

0 

0 

8.7185 

0 

a3 

0 

0 

0.4428 

0 

a4 

0 

0 

-0.1684 

0 

a5 

0 

0 

2.3099 

0 

a6 

0 

0 

-0.285 

0 

a7 

0 

0 

0.4628 

0 

a8 

0 

0 

-0.4792 

0 

a9 

0 

1 

0 

0.55 

alO 

0 

0 

0 

0.088 

all 

0 

1 

0 

0.5575 

al2 

0 

0 

0 

0.1326 

Zv( a),  %  + 

50 

55.1 

25.9 

n/a 

Zu (a),  %  - 

50 

45.5 

n/a 

23 

Total  Bound,  % 

100 

100.6 

48.9 
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Table  16.  Test  Results  For  Case  T3.C 


ai,  Naive 
Approach 

Defined 

at 

Optimized 

Max 

Optimized 

Min 

al 

0 

1 

2.4452 

0 

a2 

0 

0 

1 

0 

a3 

0 

2.5 

0.2236 

0 

a4 

0 

0 

-0.076 

0 

a5 

0 

0 

2.4456 

0 

a6 

0 

0 

2.5737 

0 

al 

0 

0 

2.5825 

0 

a8 

0 

0 

-3.1688 

0 

a9 

0 

1 

0 

1.191 

alO 

0 

0 

0 

all 

0 

1 

0 

al2 

0 

0 

0 

Zu(a),%  + 

50 

36.5 

24 

n/a 

Zv( a),  %  - 

50 

34.5 

n/a 

17.7 

Total  Bound,  % 

100 

71.0 

41.7 
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